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Authorized  for  Navy  Use  List  (NAVSEAINST  10560.2  series) 

Metric  unit  of  pressure  conveniently  sized  for  supply  pressures.  One  bar  =  100  kPa, 
or  14.5  psi. 

A  metric  expression  of  static  pressure  head.  One  cmH20  =  0.01  meters  of  pure 
water.  In  pressure  equivalents,  1  cmH20  =  0.736  torr,  981.8  Pa,  or  0.0982  kPa. 

Fraction  inspired  02.  The  fraction  of  the  inspired  gas  composed  of  oxygen. 

Feet  of  seawater,  a  unit  of  pressure.  One  fsw  =  0.3063  msw. 

Joules  per  liter,  unit  of  measure  for  "Work  of  Breathing"  normalized  for  tidal  volume. 
One  J/L  =  1  kPa. 

Kilopascals  or  newton/m2,  unit  of  pressure.  One  kPa  ~  10.2  cmH20 

Kilogram  meter  per  second2,  units  of  force,  thrust,  or  buoyancy.  From  Newton’s 
second  law:  F  =  m-a.  One  Ibf  (pound  force)  =  4.448  kg  m  sec"2  =  4.448  newton. 

Meters  of  sea  water.  One  msw  =  3.2646  fsw. 

Naval  Sea  Systems  Command 

Navy  Experimental  Diving  Unit 

Pounds  per  square  inch,  an  English  measure  of  pressure.  One  psi  =  6.895  kPa. 

1  bar  =  14.504  psi. 

Volume  averaged  pressure,  or  resistive  effort,  otherwise  known  by  the  misnomer 
Work  of  Breathing  (WOB).  A  computer-derived  estimate  of  total  resistive  respiratory 
effort  obtained  when  breathing  a  UBA  with  a  mechanical  breathing  simulator. 

Respiratory  minute  volume  with  units  of  L-min"1.  In  scientific  publications,  this  is 
referred  to  as  expired  ventilation  ( VE ). 

Standard  temperature  (0°C),  pressure  (1  atm  abs),  dry. 

Metabolic  oxygen  consumption  in  L-min'1  at  STPD  conditions. 
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INTRODUCTION 


The  U.S.  Navy  Experimental  Diving  Unit  (NEDU)  evaluates  samples  of  sodalime  CO2 
absorbent  for  use  in  diving  and  hyperbaric  applications.  The  diving  Navy  currently  uses 
both  large  and  small  grain  absorbent.  The  large  grain  absorbent,  with  granule  sizes 
ranging  from  4  to  8  U.S.  mesh,  is  the  traditional  diving  grade  absorbent.  High 
Performance  (HP)  Sodasorb,  manufactured  by  W.R.  Grace;  Sofnolime408  (Molecular 
Products,  Thaxted,  U.K.);  and  Divesorb  Pro  (Drager,  Germany)  are  such  absorbents. 
Small  grain  absorbent  granules  range  from  10  to  14  U.S.  mesh  (8-12  British  mesh  size, 
Table  1),  and  are  used  to  provide  extended  canister  durations.  Molecular  Product’s 
Sofnolime  812  is  a  typical  small  grain  absorbent. 

When  NEDU  evaluates  sodalime  quality,  we  determine,  among  other  things,  whether  or 
not  the  absorbent  meets  specifications  for  granule  size  distributions.  However,  the  U.S. 
Navy  buys  absorbent  from  at  least  three  countries  (U.S.,  U.K  and  Germany)  with  each 
supplier  having  different  sodalime  specifications.  Because  of  disparities  between  sieve 
sizes  used  in  testing  and  mesh  sizes  required  in  various  commercial  and  national 
specifications,  a  simple  comparison  is  not  always  possible. 

The  purposes  of  this  report  are  fourfold.  We  first  wish  to  understand  the  impact  of 
existing  and  proposed  specifications  on  sodalimes  with  varying  spreads  of  granule 
sizes.  Previous  reports  have  dealt  only  with  differences  in  mean  granule  size. 

Second,  we  wish  to  establish  a  rigorous  procedure  by  which  the  test  operator  can,  with 
minimal  effort,  decide  if  a  sodalime  sample  meets  specifications.  This  work  will  show 
that  sodalime  specifications  are  liberal  in  terms  of  the  allowed  distributions  of  granule 
sizes.  Consequently,  to  ensure  fairness  to  prospective  sodalime  suppliers,  a  method  is 
needed  for  comparing  sodalime  samples  other  than  just  comparing  a  new  sample  with 
existing  data. 

A  third  purpose  is  to  describe  how  a  sodalime  sample  can  be  checked  against  a  foreign 
specification  if  the  foreign  sieve  sizes  are  not  available  to  the  testing  laboratory.  We  will 
see  that  the  foreign  standards  can  be  converted  mathematically  to  a  form  any  laboratory 
can  use. 

NEDU  has  used  the  procedure  described  in  this  report  for  several  years,  and  thus  the 
fourth  purpose  of  this  report  is  to  document  the  existing  methodology.  Although  we  use 
examples  relevant  to  the  testing  sieves  available  at  NEDU,  the  procedures  are  general 
and  can  be  adapted  to  virtually  any  assemblage  of  testing  sieves  or  applied  to  any 
national  or  corporate  specification. 
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I.  DISTRIBUTION  DETERMINATION 


METHODS 

To  prevent  this  mathematical  analysis  from  becoming  too  unwieldy,  we  must  assume 
that  the  probable  granule  size  distributions  are  small  in  number.  We  hypothesized  that, 
based  on  first  principles,  the  distributions  should  be  either  normal  (Gaussian)  or  log¬ 
normal.  This  presumption  is  reasonable  due  to  the  random  nature  of  the  crushing 
process  used  to  generate  granules  from  stock  material  and  to  the  particular  applicability 
of  the  normal  and  lognormal  distributions  to  random  processes1'4. 

Although  crushed  absorbent  granules  are  screened  at  the  factory  to  produce  products 
of  varying  size  ranges,  we  anticipated  that  inefficiencies  in  the  screening  process  would 
narrow  the  granule  size  distribution  rather  than  sharply  truncate  it.  That  is,  after 
screening,  the  granule  size  distribution  would  be  likely  to  remain  smooth  rather  than 
become  discontinuous. 

Our  first  task  was  to  confirm  the  reasonableness  of  our  assumptions:  that  is,  to  measure 
the  actual  granule  size  distributions  of  various  sodalime  samples.  To  our  knowledge, 
this  has  never  been  done  with  the  degree  of  precision  required  for  this  analysis.  We 
therefore  developed  an  image  analysis  technique  using  off-the-shelf  computer  software 
to  provide  data  on  granule  size.  The  distributions  of  the  granule  sizes  were  then  found 
by  fitting  size  data  to  a  number  of  model  distributions. 

Samples  containing  a  minimum  of  800  granules  of  Sofnolime  812  and  408  (Molecular 
Products,  Thaxted,  U.K),  small  and  large  granule  size  absorbents,  respectively,  were 
spread  on  the  glass  surface  of  a  paper  photocopier  and  copied  against  a  dark 
background.  A  ruler  was  included  in  the  image  as  a  size  reference.  The  result  was  a 
high  contrast  image  with  white  granules  against  a  black  background.  The  image  was 
then  scanned  by  a  digital  scanner  (Hewlett  Packard  ScanJet  3C)  and  run  by  a  Hewlett 
Packard  Vectra  computer  using  Image  Assistant.  1.12  (Caere  Corp.)  That  image  was 
then  imported  into  a  paint  program,  and  partial  images  of  granules  on  the  periphery  of 
the  image  were  manually  deleted.  The  ruler  was  also  removed  from  the  image  after 
calibration  marks  had  been  transferred  to  the  main  body  of  the  image.  Figures  1  and  2 
show  the  result.  The  calibration  marks  in  the  upper  right  corner  of  the  figures  were  20 
millimeters  (mm)  apart. 

The  images  were  then  subjected  to  automated  image  analysis  (SigmaScan  Pro  4.0, 
Jandel  Scientific).  The  information  obtained  from  each  image  was  stored  in  a 
spreadsheet  containing  the  number  of  granules  as  well  as  the  number  of  pixels  and  the 
feret  diameter  for  each  granule.  The  feret  diameter  is  the  diameter  of  a  circle  that  yields 
the  same  number  of  pixels  as  the  actual  number  of  pixels  for  each  granule  image.  In 
other  words,  ferets  convert  irregularly  shaped  images  into  equivalent  circles  with  the 
same  area.  The  calibration  marks  were  used  to  convert  the  feret  diameters  from 
numbers  of  pixels  to  mm. 
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After  the  spreadsheet  was  imported  into  the  graphing  program  SigmaPlot  (Jandel 
Scientific),  a  mathematical  transform  was  used  to  convert  the  spreadsheet  data  into  a 
histogram  representing  the  frequency  of  occurrence  of  50  feret  diameter  categories. 
Another  Jandel  software  product,  TableCurve  2D  v.  4,  fit  up  to  74  peak  distributions  to 
the  histograms.  The  curve  fits  were  accomplished  by  the  method  of  least  squares  and 
ranked  by  the  value  of  the  F  statistic  for  the  fit. 

RESULTS 

Sofnolime  812 


The  sample  of  fine  grain  Sofnolime  812  contained  1230  granules  (Figure  1).  The 
distribution  of  feret  diameters  is  shown  in  Figure  3.  The  mean  granule  size  was  1 .76 
mm.  The  peaks  occurring  in  Figure  3  at  very  small  diameter  sizes  represent  dust  and 
granule  fragments.  Since  we  had  no  interest  in  the  overall  bimodal  distribution  (granules 
plus  dust);  we  excluded  the  dust  peaks  from  the  distribution  fitting  process. 

The  best  fit  for  the  fine  grain  Sofnolime  812  absorbent  was  to  the  Complementary  Error 
Function,  closely  related  to  the  Gaussian.  The  second  best  fit  was  a  Gaussian 
distribution  with  a  mean  (p)  of  1 .73  ±0.012  mm  (best  estimate  ±  standard  error  of  the 


Figure  1.  Sofnolime  812  Granules. 
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Figure  2.  Sofnolime  408  Granules. 


4 


80 


70  - 

60  - 


N 


50  - 

40  - 

30  - 


20  A 


10 

0 


0 


1230 


1  2  3 

feret  diameter  (mm) 


Figure  3.  Granule  size  distributions  for  Sofnolime  812. 


Figure  4  shows  the  best  fit  of  the  feret  size  frequency  distribution  to  the  Gaussian  (top 
panel);  the  middle  panel  shows  fit  residuals.  Both  95%  confidence  and  prediction  limits 
for  the  fit  to  the  Gaussian  are  in  the  bottom  panel.  Comparable  graphs  for  the 
complementary  error  function  peak  are  shown  in  Figure  5. 

Gaussian  Function 


The  Gaussian  function,  Equation  (1),  is  both  continuous  and  integrable. 


y(x) 


(1) 


The  normal  or  Gaussian  distribution  that  results  from  the  application  of  this  function  is 
given  by 


-(x-M)2 

y(x )  =  e  2c7" 


(2) 


where  x  is  granule  feret  diameter  in  mm,  p  is  the  mean  granule  diameter  in  mm,  and  a 
is  the  standard  deviation  of  the  distribution. 

When  equation  (2)  is  divided  by  crV^r  ,  we  obtain  the  so-called  standard  normal  curve 
described  in  every  statistics  book.  The  area  under  the  standard  normal  (the  integral  for 
all  x)  is  equal  to  1 . 
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Figure  4.  Gaussian  fit  to  Sofnolime  812  granule  size  distribution.  The 
fit,  residual  error,  and  95%  confidence  and  prediction  limits  are  shown 
(top  to  bottom). 
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Figure  5.  Fit  of  the  peak  of  the  complementary  error  function  to 
Sofnolime  812  granule  distribution. 
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Figure  6.  Granule  size  distribution  for  Sofnolime  408. 

Sofnolime  408 


The  sample  of  large  grain  Sofnolime  408  contained  864  granules,  not  including  dust 
particles  (Figure  2).  The  distribution  of  feret  diameters  is  shown  in  Figure  6;  the  mean 
granule  size  was  3.45  mm,  which  was  approximately  twice  as  large  as  the  1 .76  mm 
mean  for  Sofnolime  812.  The  best  fit  for  the  large  grain  Sofnolime  408  absorbent  was  a 
log-normal  distribution,  with  a  p  of  3.30  ±  0.022  mm  (best  estimate  ±  standard  error  of 
the  estimate,  SEE)  and  a  a  of  0. 1 95  ±  0.007  mm. 

The  term  log-normal  means  that  the  weights  of  the  granule  sizes  are  distributed 
normally  with  respect  to  the  logarithms  of  the  respective  granule  diameters.  The 
equation  for  the  log-normal  is  similar  to  that  for  the  Gaussian,  with  the  main  difference 
being  that  x  (granule  feret  diameter)  is  replaced  by  ln(x). 


~(ln(A-)-/y)2 

y(x)  =  e  2  ct2 
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II.  SPECIFICATION  TRANSLATION 


METHODS 

Specifications  for  sodalime  granule  size  do  not  describe  a  particular  distribution  of  sizes, 
but  rather  apply  bounds  to  the  possible  granule  size  distributions.  As  an  example, 
according  to  a  British  specification  for  8-12  U.K  mesh  range  material  (10-14  U.S.  mesh, 
Table  1 ),  no  more  than  20%  of  the  sieved  absorbent  may  be  smaller  than  a  U.K.  mesh 
size  of  12  (1 .40  mm).  The  specification  provides  a  very  general  requirement  that  a 
number  of  different  granule  size  distributions  could  meet. 


Table  1.  Mesh  size  conversions. 


U.S.  Sieve 

Tyler  Sieve 

British  (BS) 

ISO  Standard 

3.5 

3.5 

3 

5.60 

4 

4 

4 

4.75 

5 

5 

— 

4.0 

6 

6 

5 

3.35 

7 

7 

6 

2.80 

8 

8 

7 

2.38 

10 

— 

8 

2.00 

12 

10 

10 

1.70 

14 

12 

12 

1.40 

30 

28 

25 

0.600 

40 

35 

36 

0.425 

The  general  method  we  use  to  determine  whether  or  not  a  sodalime  sample  meets 
specification  requirements  is  by  first  characterizing  the  granule  size  distribution  of  the 
sample,  not  by  the  laborious  image  analysis  technique  described  in  Section  I  but  by  a 
sieving  procedure  more  suitable  to  large-scale  analyses. 

The  next  step  is  to  statistically  compare  the  obtained  distribution  with  a  model 
distribution  known  to  meet  the  specification.  However,  since  many  possible  model 
distributions  could  meet  the  specifications,  knowing  which  one  to  compare  with  the  test 
sample  is  not  straightforward.  An  incorrect  choice  of  a  model  distribution  could  unfairly 
penalize  a  particular  test  sample.  The  following  procedures  and  algorithms  are  designed 
always  to  select  the  best  possible  model  distribution  to  compare  with  a  particular  test 
sample. 

In  both  the  normal  (Gaussian)  and  log-normal  distributions,  the  entire  distribution  can  be 
described  by  only  two  parameters,  the  mean  (p)  and  the  standard  deviation  (a).  An 
additional  parameter  is  often  added  to  define  the  maximum  frequency  of  occurrence  in 
absolute  terms  (Appendix  A). 
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We  accept  the  premise  that  the  mean  of  the  granule  size  distribution  should  be  close  to 
the  presumed  median  of  the  mesh  range:  i.e.,  the  average  size  for  absorbent  in  the  4-8 
mesh  size  range  should  fall  in  the  vicinity  of  mesh  size  6.  In  the  results  section  our  data 
supports  that  premise.  Assuming  either  normal  or  log-normal  distributions,  we  generate 
a  range  of  candidate  distributions  with  various  values  of  p  and  a  that  center  on  the 
expected  median  granule  size.  A  range  of  os  are  examined  that  range  from 
extraordinarily  small  up  to  a  o  larger  than  would  be  expected  in  actual  samples. 

The  following  material,  including  the  appendixes,  describe  procedures  for  working  with 
the  Gaussian  distribution.  Analogous  methods  are  used  with  the  log-normal  distribution 
(Appendix  C). 

The  percentage  of  each  candidate  distribution  that  lies  within  particular  mesh  ranges 
can  be  found  by  integrating  Equation  (2)  over  definite  limits  representing  the  extremes 
of  the  respective  size  ranges  expressed  as  granule  diameter  in  millimeters  (mm)  rather 
than  mesh  size.  Therefore,  the  limits  of  integration  for  each  integral  represent  granule 
diameters  at  the  extreme  of  each  mesh  range  cited  in  the  specifications. 

Appendices  A-C  are  documents  generated  by  MathCad  (MathSoft,  Inc.,  Princeton,  NJ) 
detailing  the  steps  taken  to  select  the  few  model  distributions  that  meet  sodalime 
specifications  from  numerous  candidate  distributions.  For  both  large  grain  (4-8  U.S. 
mesh)  and  small  grain  (10-14  U.S.  mesh)  absorbent,  we  simulated  more  than  90 
candidate  distributions,  each  defined  by  a  separate  p  and  o.  For  the  example  of 
normally  distributed  10-14  mesh  sodalime,  we  simulated  99  distributions  and  tested 
them;  means  were  assumed  to  vary  between  1 .45  and  1 .95  mm  in  diameter,  with 
increments  of  0.05  mm.  Standard  deviations  varied  from  0.05  to  0.45  mm,  also  in 
increments  of  0.05  mm. 

Using  Table  1  for  mesh  size  to  mm  equivalences,  we  defined  four  integral  equations 
with  definite  limits  corresponding  to  mesh  ranges  called  out  in  the  specifications.  The 
first  step  in  assessing  the  candidate  distributions  was  to  evaluate  each  of  the  four 
integral  equations  for  all  candidate  distributions.  Each  integral  was  then  tested  against 
the  logical  conditions  listed  in  the  relevant  sodalime  specification.  For  example,  referring 
to  Appendix  A  for  fine  grain  sodalime,  and  the  matrix  A,  the  integral  equation  for  the 
Gaussian  between  limits  of  2.8  and  20  mm  (sizes  larger  than  #7  mesh),  most  but  not  all 
candidate  distributions  yielded  sample  fractions  less  than  1%  (p.  A2).  That  is,  less  than 
1  %  of  the  sodalime  sample  would  be  larger  than  2.0  mm  in  diameter.  Those 
distributions  that  met  the  specification  requirement  were  marked  by  a  1  in  the  matrix  of 
p  and  a,  and  those  that  did  not  meet  the  specification  were  indicated  by  a  0  (array  XA, 
binary  outcome). 

When  either  1  or  0  had  been  applied  to  each  cell  in  the  (p,  o)  array  for  each  integral,  the 
matrices  were  added  (xtot,  p.  A6).  Those  cells  with  a  sum  of  4  (1  for  each  integral)  were 
by  definition  cells  that  met  all  four  specification  requirements.  Those  cells  were  then 
further  identified  with  a  1  while  all  others  were  assigned  a  0  value  in  the  “Final”  array, 
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forming  a  triangular  pattern  of  ones  among  a  field  of  zeros.  Of  99  potential  distributions, 
only  46  met  the  complete 
specification  and  therefore  could  be 
classified  as  model  distributions 
(“Final”  array,  p.  B-6). 

NEDU  Sieves 

The  next  step  was  to 
reevaluate  the  normal  integrals 
over  limits  defining  the  sieve  sizes 
available  at  NEDU  (Fig.  8).  Six 
definite  integrals  representing  six 
size  ranges  were  defined  (3.5  mesh 
size  and  larger,  3.5  -  5  mesh,  5  -  8 
mesh,  8-10  mesh,  10-14  mesh, 
and  14-30  mesh).  Each  of  the  six 
integrals  was  applied  to  all  of  the  99 
potential  distributions  for  a  total  of 
594  integrations.  Obviously,  many 
of  those  results  were  of  no  interest 
since  they  applied  to  distributions 
that  did  not  meet  the  specification. 

Those  integrals  were  eliminated 
from  view  by  applying  a  masking 
function  or  filter  consisting  of  the 
Final  Array.  The  result  was  a  two- 
dimensional  array  of  p  and  o  with  a 
value  of  either  zero  or  the  integral 
for  each  of  the  six  mesh  ranges 
(Tables  2  and  3). 

Three  p  and  a  combinations  in  Table  2  represent  the  extremes  of  distributions 
meeting  the  specifications  for  10-14  mesh  sodalime.  Those  values  in  the  top  left,  bottom 
left,  and  mid-far  right  portions  of  the  table  are  in  bold  font  for  clarity.  They  represent 
small  granules  with  a  tight  distribution  of  granule  sizes,  larger  granules  with  an  equally 
tight  distribution,  or  mid-size  granules  with  a  wide  distribution,  respectively. 

The  non-zero  numbers  in  Table  2  tell  us  that  if  absorbents  meeting  the 
specifications  are  sieved,  the  #14  mesh  sieve  would  collect  anywhere  from  59%  to 
100%  of  the  test  sample,  depending  on  which  of  the  46  model  distributions  best  applied 
to  the  sample.  Both  the  shapes  of  these  three  distributions  and  their  expected  sieve 
retentions  are  shown  in  Figure  9. 
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Table  3  shows  similar  results  for  large  grain,  4-8  U.S.  mesh  sodalime.  Anywhere  from 
61%  to  79%  of  the  absorbent  could  be  trapped  on  the  #8  U.S.  mesh  sieve  and  still  meet 
specifications. 


Table  2.  The  “Final”  array  for  10-14  mesh  sodalime,  #14  sieve. 


Distribution  Standard  Deviation  (a,  mm) 

p  (mm) 

0.05 

0.10 

0.15 

0.20 

0.25 

0.30 

0.35 

0.40 

0.45 

1.45 

0.84 

0 

0 

0 

0 

0 

0 

0 

0 

1.50 

0.98 

0.84 

0 

0 

0 

0 

0 

0 

0 

1.55 

1 

0.93 

0.84 

0 

0 

0 

0 

0 

0 

1.60 

1 

0.98 

0.91 

0.83 

0 

0 

0 

0 

0 

1.65 

1 

0.99 

0.95 

0.87 

0.79 

0 

0 

0 

0 

1.70 

1 

1 

0.97 

0.9 

0.81 

0.73 

0.65 

0 

0 

1.75 

1 

1 

0.97 

0.9 

0.81 

0.73 

0.65 

0.59 

0 

1.80 

1 

1 

0.95 

0.88 

0.8 

0.72 

0.64 

0 

0 

1.85 

1 

0.98 

0.92 

0.84 

0.76 

0 

0 

0 

0 

1.90 

1 

0.95 

0.86 

0 

0 

0 

0 

0 

0 

1.95 

0.99 

0 

0 

0 

0 

0 

0 

0 

0 

Table  3.  The  “Final”  array  for  4-8  mesh  sodalime,  #8  sieve. 


Distribution  Standard  Deviation  (or,  mm) 

p  (mm) 

0.65 

0.70 

0.75 

0.80 

0.85 

0.90 

0.95 

1.0 

3.02 

0 

0 

0 

0 

0 

0 

0 

0 

3.08 

0.79 

0 

0 

0 

0 

0 

0 

0 

3.15 

0.79 

0.76 

0.73 

0 

0 

0 

0 

0 

3.21 

0.79 

0.76 

0.73 

0.69 

0 

0 

0 

0 

3.28 

0.79 

0.75 

0.72 

0.69 

0.66 

0 

0 

0 

3.34 

0.78 

0.75 

0.71 

0.69 

0.66 

0.63 

0.61 

0 

3.41 

0.77 

0.73 

0.70 

0.68 

0.66 

0.62 

0 

0 

3.47 

0.75 

0.72 

0.69 

0.66 

0.64 

0 

0 

0 

3.54 

0.73 

0.70 

0.67 

0.63 

0 

0 

0 

0 

3.60 

0.70 

0.68 

0.65 

0 

0 

0 

0 

0 

3.67 

0.67 

0.65 

0 

0 

0 

0 

0 

0 

3.73 

0.64 

0 

0 

0 

0 

0 

0 

0 
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Frequency  of  Occurence 


Sieve  Size  (mm) 


Figure  9.  Three  granule  size  distributions  (bottom  panel)  and  corresponding  sieve 
results  (top  panel). 


MeshFit 


It  is  possible  that  a  test  sample  will  match  none  of  the  model  distributions  particularly 
well.  Fortunately,  the  degree  of  match  or  mismatch  can  be  measured  and  used  as  a 
basis  for  evaluating  the  test  sample.  NEDU  uses  the  MeshFit  software  developed  in- 
house  to  enter  the  sieve  data  from  the  test  sample  and  to  locate  the  one  model 
distribution  that  best  agrees  with  the  test  sample  data.  The  goodness  of  fit  between  the 
test  data  and  the  best  model  distribution  is  then  found  through  two  mechanisms,  one 
being  simply  the  sum  of  the  absolute  differences  between  actual  and  ideal  retained 
fractions  and  the  other  being  the  classical  chi  square  statistical  test  for  goodness  of  fit. 
MeshFit  was  written  in  Visual  Basic  Version  6  for  Windows  98,  32-bit  operating  system. 
The  most  relevant  portions  of  code  are  provided  in  Appendix  D. 

As  explained  in  Appendix  G,  the  user  selects  whether  the  test  sample  is  from  a  4-8  U.S. 
mesh  sample  or  10-14  U.S.  mesh  sample.  The  user  then  enters  the  mean  retained 
percentages  for  all  sieve  screens  obtained  from  the  sieving  of  five  or  more  sodalime 
samples.  The  software  next  identifies  the  best  fit  to  a  model  distribution  and  plots  bar 
graphs  for  the  mesh  distribution  of  both  the  sample  and  the  best  model.  Finally,  the  chi 
square  goodness  of  fit  statistic  is  computed,  with  its  significance  being  determined  from 
the  chi  square  and  the  degrees  of  freedom  for  the  test. 

Chi  square  Goodness  of  Fit 

The  chi  square  test  enables  observed  versus  expected  values  for  multiple  data 
groupings  to  be  tested.  In  this  particular  case,  we  compare  the  observed  fraction 
retained  by  a  particular  sieve  screen  against  the  expected  fraction  for  the  closest 
distribution,  either  normal  or  log-normal,  meeting  the  specifications. 

Although  the  expected  values  may  be  real  numbers,  to  be  valid  the  chi  square  test 
requires  that  all  observations  be  expressed  as  integers.  Furthermore,  to  be  included  in 
the  test,  the  minimum  value  observed  must  be  5  or  greater.  Data  groupings;  i.e.,  mesh 
sizes,  may  be  pooled  to  meet  the  requirement  of  5  observations.  To  satisfy  the  above 
requirements,  when  the  chi  square  test  is  run  retained  fractions  are  expressed  as 
percentages  of  the  total  sample. 


RESULTS 

MeshFit 

Figures  10  and  1 1  are  MeshFit  program  printouts  illustrating  the  analysis  of  test  results 
from  two  large  grain  sodalime  samples.  In  Figure  10  the  sample  mesh  distribution  was 
statistically  equivalent  to  the  best  model  distribution  for  4-8  mesh  absorbent.  In  Figure 
1 1  the  sample  distribution  departed  significantly  from  the  expected  mesh  distribution. 
The  probability  that  a  sample  distribution  would  depart  as  far  as  it  did  from  the  model  if 
in  fact  the  sodalime  met  the  best  model  distribution  was  <  0.001 .  Statistical  significance 
is  assumed  at  a  probability  (P)  less  than  0.05. 
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Figure  10.  Screen  print  from  MeshFit,  4-8  mesh  sodalime  data.  The  measured  and  model  distribution  characteristics  are 
similar. 
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Figure  11.  MeshFit  screen  print;  4-8  mesh  data  not  matching  any  of  the  allowed  distributions.  The  closest  match  is 
displayed. 
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Appendix  D  results  from  printing  “details”  when  running  the  MeshFit  program  for  the 
sample  in  Figure  10.  The  appendix  lists  the  results  of  fitting  the  sample  distribution  to 
each  of  the  96  candidate  distributions.  The  first  column  records  the  number  of  the 
candidate  distribution,  the  second  presents  the  error  of  the  fit,  and  the  remaining 
columns  show  the  expected  fractions  of  retained  sodalime  on  the  3.5-,  5-,  8-,  10-,  14-, 
and  30-mesh  screens.  Only  the  distributions  that  have  a  non-zero  value  in  two  or  more 
of  the  retained  fraction  columns  serve  as  model  distributions.  As  shown  in  Figure  10, 
the  best  fit  is  obtained  with  model  distribution  number  59.  The  expected  retained 
fractions  in  the  Details  printout  for  row  59  (marked  in  Appendix  D  by  a  shaded  line)  are 
also  displayed  in  Figure  10  under  the  heading  “Model  Mesh  Fractions”.  As  expected, 
the  mean  for  the  best  model  distribution  corresponds  to  a  mesh  size  of  5.9, 
approximately  midway  between  the  mesh  sizes  of  4  and  8.  The  majority  of  the  granules 
have  a  mesh  size  of  5.1  (mean  -  SD)  -  7.2  (mean  +  SD). 

At  the  end  of  the  “Details”  listing  is  an  across  the  board  comparison  between  the 
observed  results  (obs)  and  the  model  distribution  (modi).  Retained  fractions  from  both 
sample  and  model  populations  have  been  pooled  as  necessary  to  meet  the 
requirements  of  the  chi  square  test:  to  have  a  minimum  of  “5”  in  any  of  the  compared 
groupings. 

In  the  example  of  Figure  10,  only  three  pooled  mesh  sizes  (n)  can  be  used  for  the  chi 
square  test.  The  resulting  degrees  of  freedom  for  the  test  are  n-1 ,  or  2,  as  reported  in 
Figure  10.  The  probability  of  obtaining  a  chi  square  goodness  of  fit  statistic  of  0.01 63 
from  a  distribution  that  is  the  same  as  the  model  distribution  is  0.998.  Since  that 
probability  is  larger  than  0.05,  we  accept  the  null  hypothesis  that  the  tested  sample  and 
the  model  are  equivalent.  That  information  is  displayed  on  the  MeshFit  screen  by  the 
wording  "Same  as  model",  or  in  an  alternate  version  of  the  software  (Appendix  G), 
“Meets  spec.” 

Appendix  E  is  the  “Detail”  listing  for  Figure  1 1 .  The  test  sample  fit  candidate  distribution 
89  (row  89,  shaded)  the  best.  This  time  after  pooling,  only  2  mesh  sizes  could  be  tested 
by  the  chi  square;  therefore,  the  test  had  1  degree  of  freedom.  The  probability  of 
obtaining  a  chi  square  statistic  of  20.65  from  a  test  sample  that  matches  the  model  is 
less  than  0.001 .  Therefore,  we  reject  the  null  hypothesis  that  the  test  sample  matches 
the  best  modeled  distribution;  that  is,  the  sample  differs  significantly  from  the  model  and 
therefore  does  not  meet  specifications.  That  information  is  displayed  on  the  MeshFit 
screen  by  the  wording  "Significantly  different  from  best  model!"  or  in  Appendix  G,  “Does 
not  meet  spec!” 


DISCUSSION 

Image  analysis  of  absorbent  granules  show  that  the  assumptions  required  for  the 
MeshFit  program  are  reasonable.  Both  normal  and  log-normal  distributions  described 
the  absorbent  granule  size  distributions  of  our  two  samples  (Figures  1-7).  Furthermore, 
the  mean  granule  size  did  lie  near  the  middle  of  the  mesh  range  of  the  absorbent. 
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The  NEDU  image  analysis  technique  for  determining  granule  size  distributions  is  too 
time  consuming  for  routine  use.  The  standard  sieve  test  with  Rotap  shaker  is  far  more 
efficient.  However,  armed  with  the  insights  gained  from  image  analysis,  we  can 
mathematically  analyze  sieve  data  to  achieve  our  goal  of  establishing  a  rigorous 
procedure  for  qualifying  or  disqualifying  a  sodalime  sample  while  requiring  minimal 
decisions  by  the  test  operator.  The  means  to  that  end  is  NEDU's  MeshFit  software, 
which  requires  as  input  both  the  mean  data  from  multiple  sieving  runs  and  the  data  files 
generated  by  MathCad  (Appendices  A-C). 

One  benefit  from  this  approach  is  that  this  method  reduces  unfairness  in  comparing  a 
sodalime  sample  to  specifications.  The  approach  appreciates  that  specifications  are 
fairly  general.  If  MeshFit  says  that  a  sodalime  sample  does  not  meet  specification,  it  is 
satisfying  to  know  that  the  failing  sample  has  been  compared  to  virtually  all  possible 
distributions  meeting  the  specification.  Mesh  Fit’s  accept  or  reject  decision  is  based  on  a 
comparison  with  the  test  sample's  best  match. 

One  of  the  aims  of  this  work  was  to  learn  how  existing  and  proposed  specifications  for 
sodalimes  would  be  affected  by  varying  spreads  of  granule  sizes.  As  shown  in  Figure  9, 
current  specifications  for  fine  grain  absorbent  allow  widely  varying  distributions  of 
granule  sizes.  Narrow  distributions  can  differ  considerably  in  their  mean  value  and  still 
meet  the  specification.  Since  mean  granule  size  is  inversely  related  to  absorbent 
activity,  which  in  turn  dramatically  affects  canister  duration,  specifications  for  absorbents 
with  narrowly  defined  granule  sizes  should  have  much  tighter  restrictions  on  allowable 
size  than  specifications  for  absorbents  with  widely  dispersed  granule  sizes. 
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APPENDIX  A.  Normally  Distributed  10-14  Mesh  Sodalime 

Problem:  A)  Find  a  Gaussian  (normal)  distribution  that  fulfills  the  requirements  of  the  original 
specification  for  8-12  (10-14  US)  mesh  size  C02  absorbent  (sodalime). 

B)  Find  allowable  %  retention  of  absorbent  for  sieve  sizes  used  in  the  NEDU  T&E  laboratory. 


ORIGIN*  1 

j  :=  1 ..  1 1  i  :=  1 ..  9  Iteration  range 

Gaussian  (normal)  distribution  parameters 
Mean  Standard  deviation 

Mj:=  1.40+ .05j  oj^i-.OS 


Normal  distribution: 


y(x,p  ,o):= 


(»-u)2 

[2-(o)2] 


Integral  equations  relevant  to  specifications  for  10-14  mesh  absorbent 
All  mesh  sizes  expressed  in  diameter  of  mesh  screen  (in  mm). 
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The  proposed  draft  sodalime  specification  calls  for  no  more  than  1%  larger  than  7  mesh,  no  more  than 
30%  in  the  7-10  mesh  size  range,  no  less  than  48%  in  the  10-14  range,  no  more  than  20%  between  14 
and  40  mesh. 

Out  of  99  candidate  distributions,  we  found  46  that  met  all  criteria  for  the  specification,  and  were 
therefore  classified  as  model  distributions.  The  manner  in  which  these  46  distributions  were  identified  are 
indicated  on  the  following  pages.  Of  those  46,  6  are  listed  on  the  last  two  pages. 


A-1 


Fraction  of  absorbent  larger  than  #7  mesh 


Matrix  of  absorbent  fraction  greater  than 
#7  mesh  (2.8  mm)  in  size  for  various  distribution 
j's  (means)  and  i's  (standard  deviations). 
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Specific  Specification  Criterion 
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1  =  meets  specification 
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Each  ”1"  identifies  a  theoretical  distribution  that 
meets  the  portion  of  the  specification  that 
requires  no  more  than  1%  of  the  absorbent  can 
be  larger  than  #7  mesh. 
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Fraction  of  absorbent  between  #7  and  #10  mesh 
(2.8  mm  to  2.0  mm) 
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l 

1 

1 

1 

I 

1 

1 

1 

1 

i 

1 

1 

1 

1 

1 

1 

1 

1 

l 

1 

1 

1 

1 

1 

1 

1 

1 

l 

1 

1 

1 

1 

1 

1 

1 

1 

l 

1 

1 

1 

1 

1 

1 

1 

1 

l 

1 

1 

1 

1 

0 

0 

1 

1 

l 

1 

1 

0 

0 

0 

0 

1 

1 

l 

0 

0 

0 

0 

0 

0 

U 

0 

0 

0 

0 

0 

0 

0 

0, 

Specific  Specification  Criterion 
XB.  j:=  if(B.  .<-30,1,0) 


A- 3 


Fraction  of  absorbent  between  #10  and  #14  mesh 
(2.0  mm  to  1.4  mm) 


© 

bo 

0.69 

0.63 

0.6 

0.57 

0.53 

0.5 

0.47 

0.43 

0.98 

0.84 

0.75 

0.69 

0.63 

0.58 

0.54 

0.49 

0.45 

1 

0.93 

0.84 

0.76 

0.69 

0.62 

0.57 

0.52 

0.47 

1 

0.98 

0.9 

0.82 

0.73 

0.66 

0.59 

0.53 

0.48 

1 

0.99 

0.94 

0.85 

0.76 

0.68 

0.6 

0.54 

0.49 

1 

1 

0.95 

0.87 

0.77 

0.68 

0.61 

0.55 

0.5 

1 

0.99 

0.94 

0.85 

0.76 

0.68 

0.6 

0.54 

0.49 

1 

0.98 

0.9 

0.82 

0.73 

0.66 

0.59 

0.53 

0.48 

1 

0.93 

0.84 

0.76 

0.69 

0.62 

0.57 

0.52 

0.47 

0.98 

0.84 

0.75 

0.69 

0.63 

0.58 

0.54 

0.49 

0.45 

,0.84 

0.69 

0.63 

0.6 

0.57 

0.53 

0.5 

0.47 

0.43 

Binary  Outcome 


f  j 

1 

1 

1 

1 

1 

1 

0 

0N 

1 

1 

1 

1 

1 

1 

1 

1 

0 

1 

1 

1 

1 

1 

1 

1 

1 

0 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

I 

1 

1 

1 

1 

0 

1 

1 

1 

1 

1 

1 

1 

1 

0 

J 

1 

1 

1 

1 

1 

1 

0 

o. 

Specific  Specification  Criterion 

XC.  .:=  iffC.  .>.48,1, O’) 
j.i  \  J.'  ) 


A-4 


Fraction  of  absorbent  between  #14  and  #30  mesh 
(1.40  mm  to  0.60  mm) 


0.16 

0.31 

0.37 

0.4 

0.42 

0.43 

0.44 

0.45 

0.44  N 

0.02 

0.16 

0.25 

0.31 

0.34 

0.37 

0.39 

0.4 

0.4 

1.35  x  10“ 

3 

0.07 

0.16 

0.23 

0.27 

0.31 

0.33 

0.35 

0.36 

3.36  x  10' 

5 

0.02 

0.09 

0.16 

0.21 

0.25 

0.28 

0.31 

0.32 

3.28  x  10" 

7 

6.21  x  10" 

3 

0.05 

0.11 

0.16 

0.2 

0.24 

0.26 

0.29 

1.24  x  10" 

9 

1.35  x  10" 

3 

0.02 

0.07 

0.12 

0.16 

0.2 

0.23 

0.25 

1.79  x  10" 

12 

2.33  x  10" 

4 

9.82  x  10" 

3 

0.04 

0.08 

0.12 

0.16 

0.19 

0.22 

0 

3.17  x  10" 

5 

3.83  x  10" 

3 

0.02 

0.05 

0.09 

0.13 

0.16 

0.19 

0 

3.41  x  10" 

6 

1.35  x  10" 

3 

0.01 

0.04 

0.07 

0.1 

0.13 

0.16 

0 

2.88  x  10" 

7 

4.29  x  10" 

4 

6.21  x  10'  3 

0.02 

0.05 

0.08 

0.11 

0.13 

0 

1.92  x  10" 

8 

1.23  x  10“ 

4 

2.98  x  10"  3 

0.01 

0.03 

0.06 

0.08 

0.11  , 

Binary  Outcome 

Specific  Specification  Criterion 
XD., 


/j 

0 

0 

0 

0 

0 

0 

0 

0N 

1 

1 

0 

0 

0 

0 

0 

0 

0 

1 

1 

1 

0 

0 

0 

0 

0 

0 

1 

1 

1 

1 

0 

0 

0 

0 

0 

1 

1 

1 

1 

1 

0 

0 

0 

0 

1 

1 

1 

1 

1 

1 

1 

0 

0 

1 

1 

1 

1 

1 

1 

1 

1 

0 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

\1 

1 

1 

1 

1 

1 

1 

1 

1 J 

A-5 


(A4) 


r  4 

3 

3 

3 

3 

3 

3 

2 

2  > 

4 

4 

3 

3 

3 

3 

3 

3 

2 

4 

4 

4 

3 

3 

3 

3 

3 

2 

4 

4 

4 

4 

3 

3 

3 

3 

3 

4 

4 

4 

4 

4 

3 

3 

3 

3 

XTOT  = 

4 

4 

4 

4 

4 

4 

4 

3 

3 

4 

4 

4 

4 

4 

4 

4 

4 

3 

4 

4 

4 

4 

4 

4 

4 

3 

2 

4 

4 

4 

4 

4 

3 

3 

3 

1 

4 

4 

4 

3 

3 

3 

3 

2 

1 

,4 

3 

3 

3 

3 

3 

3 

1 

1, 

Sum  of  binary  outcomes 

XTOT.  .  :=  fXA,  .  +  XB.  .  +  XC.  .  +  XD. 
J.i  l  J.i  J.'  J.i  1 


FINAL.  .  :=  if  (XTOT.  .  =  4 , 1 ,0) 

1's  identify  distributions  that  simultaneously 
meet  all  four  specification  requirements. 
That  is,  these  are  good  model  distributions. 

This  matrix  is  used  below  as  a  mask  or 
filter  to  select  valid  values  for  retained 
fractions  with  NEDU’s  mesh  screens. 


'1 

1 


FINAL  = 


0000000  0^ 
1  0  0  0  0  0  0  0 

1  1  0  0  0  0  0  0 

1  1  1  0  0  0  0  0 

1  1  1  1  0  0  0  0 

1111110  0 
11111110 
1111110  0 
11  1  1  0  0  0  0 
1  1  0  0  0  0  0  0 

0000000  0, 


A-6 


The  following  definite  integrals  define  the  fraction  of  the  normal  distribution  falling  within 
the  confines  of  the  six  sieves  available  at  NEDU. 

3.5  sieve 

5  sieve 

8  sieve 

E.  .:= 

J 

r20 

dl 

5.6 

z-5.6 

Fi  i :=  y(x,Hj,aj)dx 

J’  J4.0 

r4.0 

G:  i:=  y(*.Mj.Oi)  d* 

J2.36 

10  sieve 

14  sieve 

30  sieve 

H.  .:= 
j»*  j 

/•2.36 

y(*.f»j.<*i)  d* 

2.06 

r2.06 

K  :=  y(x,  Hj  .CTj)  dx 

J>  J1.4 

fl.4 

L.  ;:=  y(*.Mj.Oj)dx 

*  .6 

A-7 


3.5  sieve 


E  = 


''o  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 


0  0  0  0  0 
0  0  0  0  0 
0  0  0  0  0 
0  0  0  0  0 
0  0  0  0  0 
0  0  0  0  0 
0  0  0  0  0 
0  0  0  0  0 
0  0  0  0  0 
0  0  0  0  0 
0  0  0  0  0 


0N 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0, 


Fractions  retained  by  a  #3.5  mesh 
sieve. 


Results  of  binary  filtering 


EANS  = 


'o  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
,0  0  0 


0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 


0  0  0^ 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  Oj 


Apply  binary  filter 

EANS.  .:=  E.  -FINAL.  . 
J.r  J,l  J.i 


Write  results  to  a  data  file  for 
use  by  MeshFit  software. 

VVRITEPRN(”meshl014.pm" )  :=  EANS 


A-8 


5  sieve 


( 

0 

0 

0 

0 

0 

0 

1.6  x  10“ 

13 

9.15  x  10“  11 

7.28  x  10“  9 

0 

0 

0 

0 

0 

0 

4.57  x 

10" 

13 

2.05  x  10"  10 

1.38  x  10"8 

0 

0 

0 

0 

0 

0 

X 

oo 

10" 

12 

4.53  x  10"  10 

2.6  x  10“  8 

0 

0 

0 

0 

0 

0 

3.51  x  10 

12 

9.87  x  10"  10 

4.82  x  10"8 

0 

0 

0 

0 

0 

2.38  X  10"  15 

9.45  x 

10“ 

12 

2.11  x  10"9 

8.84  x  10"8 

F  = 

0 

0 

0 

0 

0 

8.85  x  10"  15 

2.49  x  10“ 

11 

4.46  x  10"  9 

1.6  x  10"  7 

0 

0 

0 

0 

0 

3.2  x  10“  14 

6.44  x 

10“ 

11 

9.28  x  10"  9 

2.87  x  10"  7 

0 

0 

0 

0 

0 

1.12  x  10"  13 

1.63  x 

10" 

10 

1.9  x  I0"8 

5.07  x  10"7 

0 

0 

0 

0 

0 

3.85  x  10"  13 

4.05  x 

10“ 

10 

3.83  x  10"  8 

8.86  x  10"7 

0 

0 

0 

0 

0 

1.28  x  10"  12 

9.87  x 

10“ 

10 

7.61  x  10"8 

1.53  x  10“  6 

kO 

0 

0 

0 

0 

4.15  x  10"  12 

2.35  x 

10“ 

9 

1.49  x  10"  7 

2.61  x  10"  6 

Results  of  Binary  filtering 


Apply  binary  filter 


^0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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0 

0 

0 

FANS  = 

0 

0 

0 

0 

0 

8.85  X  10“  15 

2.49  x  10' 

11 

0 

0 

0 

0 

0 

0 

0 

-  14 

3.2x10 

6.44  x  10" 

11 

9.28  x  10"9 

0 

0 

0 

0 

0 

0 

1.12  x  10"13 

1.63  x  10" 

10 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

,0 

0 

0 

0 

0 

0 

0 

0 

0 

FANS.  .:=  F.  .  FINAL.  . 
j.i  i.l  J.i 

Append  results  to  the  data  file. 
APPENDPRN("meshl014.pm")  :=  FANS 

Maximum  integral 
max(FANS)  =  9.28  x  l<f  9 


A-9 


8  sieve 


/ 

0 

0 

6.73  x 

O 

1 

o 

2.69  x  10"  6 

1.36  x  10" 

4 

1.21  x  10~3 

1  4.66  x  10-3 

0.01 

0.02 

0 

0 

5.04  x 

10- 9 

8.55  x  10"  6 

2.91  x  10" 

4 

2.07  x  10"3 

1  7x  10"  3 

0.02 

0.03 

0 

0 

3.39  x 

10- 8 

2.56  x  10-5 

5.98  x  10“ 

4 

3.47  x  10~3 

0.01 

0.02 

0.04 

0 

-  14 

1.94  x  10 

2.05  x 

10"  7 

7.24  x  10"  5 

1.18  x  10" 

3 

5.65  x  10-3 

0.01 

0.03 

0.05 

0 

7.83  x  10“  13 

1.11  x 

10-6 

1.93  x  10"  4 

2.26  x  10" 

3 

8.97  x  10"  3 

0.02 

0.04 

0.06 

G  = 

0 

2.48  x  10"  " 

5.44  x 

10- 6 

4.83  x  10"  4 

4.15  x  10" 

3 

0.01 

0.03 

0.05 

0.07 

0 

6.14  x  10"  10 

2.39  x 

10-5 

1.14  x  10-3 

7.34  x  10" 

3 

0.02 

0.04 

0.06 

0.09 

0 

1.2  x  10'8 

9.46  x 

10-5 

2.55  x  10"  3 

0.01 

0.03 

0.05 

0.08 

0.11 

0 

1.84  x  10-7 

3.37  x 

io" 4 

5.39  x  10-3 

0.02 

0.04 

0.07 

0.1 

0.13 

0 

2.23  x  10-  6 

1.08  x 

10- 3 

0.01 

0.03 

0.06 

0.09 

0.13 

0.15 

<0 

2.13  x  10“  5 

3.13  x 

10- 3 

0.02 

0.05 

0.09 

0.12 

0.15 

0.18 

'0 

0 

0 

0 

0 

0 

0 

0 

0^ 

Apply  binary  filter 

0 

0 

0 

0 

0 

0 

0 

0 

0 

GANS.  .  :=  G.  .  FINAL. 
J.'  J.'  J 

0 

0 

3.39  x  10~8 

0 

0 

0 

0 

0 

0 

0 

-  14 

1.94  x  10 

2.05  x  10"  7 

7.24  x  10"  5 

0 

0 

0 

0 

0 

0 

7.83  x  10“  13 

1.11  x  IO"6 

1.93  x  10"  4 

2.26  x  IO"3 

0 

0 

0 

0 

CANS  = 

0 

2.48  x  10"  11 

5.44  x  IO"6 

4.83  x  IO-4 

4.15  x  IO-3 

0.01 

0.03 

0 

0 

0 

6.14  x  10"  10 

2.39  x  10"  5 

1.14  x  10"  3 

7.34  x  10"  3 

0.02 

0.04 

0.06 

0 

Maximum  integral 

0 

1.2  x  10"  8 

9.46  x  IO"5 

2.55  x  10"  3 

0.01 

0.03 

0.05 

0 

0 

max  (GANS)  =  0.06 

0 

1.84  x  10"  7 

3.37  x  IO"4 

5.39  x  10"  3 

0.02 

0 

0 

0 

0 

0 

2.23  x  10“  6 

1.08  x  10"  3 

0 

0 

0 

0 

0 

0 

vO 

0 

0 

0 

0 

0 

0 

0 

Oj 

A-10 


10  sieve 


0 

5.3  x  10"  10 

2.38  x  10” 

5 

1 .14  x  10"3 

7.21  x  10"  3 

0.02 

0.04 

0.05 

N 

0.07 

0 

1.07  x  10"  8 

9.45  x  10” 

5 

2.55  x  10”  3 

0.01 

0.03 

0.05 

0.06 

0.08 

0 

1.7  x  10"7 

3.37  x  10” 

4 

5.36  x  10”  3 

0.02 

0.04 

0.06 

0.08 

0.09 

0 

2.11  x  10“  6 

1.08  x  10” 

3 

0.01 

0.03 

0.06 

0.08 

0.1 

0.11 

0 

2.07  x  10"  5 

3.13  x  10” 

3 

0.02 

0.05 

0.08 

0.1 

0.11 

0.12 

H  = 

3.02  x  10"  13 

1.59  x  10"4 

8.19  x  10” 

3 

0.04 

0.07 

0.1 

0.12 

0.13 

0.14 

2.83  x  10”  10 

9.68  x  10”  4 

0.02 

0.06 

0.1 

0.13 

0.15 

0.16 

0.16 

9.97  x  10"  8 

4.66  x  10"  3 

0.04 

0.09 

0.14 

0.16 

0.17 

0.18 

0.18 

1.33  x  lO-5 

0.02 

0.08 

0.14 

0.18 

0.2 

0.2 

0.2 

0.19 

6.87  x  10"  4 

0.05 

0.14 

0.2 

0.23 

0.23 

0.23 

0.22 

0.21 

,  0.01 

0.14 

0.23 

0.27 

0.28 

0.27 

0.26 

0.24 

0.22, 

HANS  = 


0 

0 

0 

0 

0 

0 

0 

0 

O'1 

0 

1.07  x  10”8 

0 

0 

0 

0 

0 

0 

0 

0 

1.7  x  10"  7 

3.37  x  10”  4 

0 

0 

0 

0 

0 

0 

0 

2.11  x  10”6 

1.08  x  10”3 

0.01 

0 

0 

0 

0 

0 

0 

2.07  x  10”5 

3.13  x  10”3 

0.02 

0.05 

0 

0 

0 

0 

3.02  x  10”  13 

1.59  x  10”  4 

8.19  x  10”3 

0.04 

0.07 

0.1 

0.12 

0 

0 

2.83  x  10”  10 

9.68  x  I0”4 

0.02 

0.06 

0.1 

0.13 

0.15 

0.16 

0 

9.97  x  10”8 

4.66  x  10”  3 

0.04 

0.09 

0.14 

0.16 

0.17 

0 

0 

1.33  x  10”5 

0.02 

0.08 

0.14 

0.18 

0 

0 

0 

0 

6.87  x  10”  4 

0.05 

0.14 

0 

0 

0 

0 

0 

0 

0.01 

0 

0 

0 

0 

0 

0 

0 

o. 

Apply  binary  filter 

HANS.  .:=  H.  .  FINAL, 
i.i  j>  ■  J. 


Maximum  integral 

max  (HANS)  =  0.18 


A-1 1 


14  sieve 


'0.84 

0.69 

0.63 

0.6 

0.57 

0.55 

0.52 

0.49 

0.46 

0.98 

0.84 

0.75 

0.69 

0.64 

0.6 

0.56 

0.52 

0.48 

1 

0.93 

0.84 

0.77 

0.71 

0.65 

0.59 

0.55 

0.5 

1 

0.98 

0.91 

0.83 

0.76 

0.68 

0.62 

0.57 

0.52 

1 

0.99 

0.95 

0.87 

0.79 

0.71 

0.64 

0.58 

0.53 

1 

1 

0.97 

0.9 

0.81 

0.73 

0.65 

0.59 

0.54 

1 

1 

0.97 

0.9 

0.81 

0.73 

0.65 

0.59 

0.54 

1 

1 

0.95 

0.88 

0.8 

0.72 

0.64 

0.58 

0.53 

1 

0.98 

0.92 

0.84 

0.76 

0.69 

0.63 

0.57 

0.52 

1 

0.95 

0.86 

0.78 

0.72 

0.66 

0.6 

0.55 

0.51 

v0.99 

0.86 

0.77 

0.71 

0.66 

0.61 

0.57 

0.52 

0.49 

'0.84 

0 

0 

0 

0 

0 

0 

0 

o'' 

0.98 

0.84 

0 

0 

0 

0 

0 

0 

0 

Apply  binary  filter 

1 

0.93 

0.84 

0 

0 

0 

0 

0 

0 

KANS.  .:=  K.  .  FINAL 

JJ  J.> 

1 

0.98 

0.91 

0.83 

0 

0 

0 

0 

0 

1 

0.99 

0.95 

0.87 

0.79 

0 

0 

0 

0 

KANS  = 

1 

1 

0.97 

0.9 

0.81 

0.73 

0.65 

0 

0 

1 

1 

0.97 

0.9 

0.81 

0.73 

0.65 

0.59 

0 

1 

1 

0.95 

0.88 

0.8 

0.72 

0.64 

0 

0 

1 

0.98 

0.92 

0.84 

0.76 

0 

0 

0 

0 

Maximum  integral 

0.95 

0.86 

0 

0 

0 

0 

0 

0 

max  (KANS)  =  1 

,0.99 

0 

0 

0 

0 

0 

0 

0 

0, 
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30  sieve 


0.16 

0.31 

0.37 

0.02 

0.16 

0.25 

1.35  x  10-3 

0.07 

0.16 

3.24  x  10' 5 

0.02 

0.09 

3.06  x  10“  7 

6.21  x  10"  3 

0.05 

1.12  x  10”  9 

1.35  x  10“  3 

0.02 

1.57  x  10“  12 

2.33  x  10“  4 

9.82  x  10“ 

0 

3.17  x  10“  5 

3.83  x  10“ 

0 

3.4  x  10“6 

1.35  x  10“ 

0 

2.87  x  10“  7 

4.29  x  10“ 

0 

1.91  x  10“  8 

1.23  x  10“ 

0.4 

0.42 

0.43 

0.44 

0.43 

0.43  > 

0.31 

0.34 

0.37 

0.38 

0.39 

0.39 

0.23 

0.27 

0.31 

0.33 

0.35 

0.35 

0.16 

0.21 

0.25 

0.28 

0.3 

0.32 

0.11 

0.16 

0.2 

0.24 

0.26 

0.28 

0.07 

0.12 

0.16 

0.19 

0.22 

0.25 

0.04 

0.08 

0.12 

0.16 

0.19 

0.21 

0.02 

0.05 

0.09 

0.13 

0.16 

0.18 

0.01 

0.04 

0.07 

0.1 

0.13 

0.16 

6.21  x  10“  3 

0.02 

0.05 

0.08 

0.11 

0.13 

2.98  x  10' 3 

0.01 

0.03 

0.06 

0.08 

0.11  , 

f  0.16 

0 

0 

0 

0 

0 

0 

0 

o' 

0.02 

0.16 

0 

0 

0 

0 

0 

0 

0 

1.35  x  10“  3 

0.07 

0.16 

0 

0 

0 

0 

0 

0 

3.24  x  10“  5 

0.02 

0.09 

0.16 

0 

0 

0 

0 

0 

3.06  x  10“  7 

6.21  x  10' 3 

0.05 

0.11 

0.16 

0 

0 

0 

0 

LANS  = 

1.12  x  10' 9 

1.35  x  10' 3 

0.02 

0.07 

0.12 

0.16 

0.19 

0 

0 

1.57  x  10“  12 

2.33  x  10"  4 

9.82  x  10“  3 

0.04 

0.08 

0.12 

0.16 

0.19 

0 

0 

3.17  x  10' 5 

3.83  x  10“  3 

0.02 

0.05 

0.09 

0.13 

0 

0 

0 

3.4  x  10'6 

1.35  x  10“  3 

0.01 

0.04 

0 

0 

0 

0 

0 

2.87  x  10“  7 

4.29  x  10“  4 

0 

0 

0 

0 

0 

0 

,  o 

0 

0 

0 

0 

0 

0 

0 

0, 

Apply  binary  filter 

LANS.  .  :=  L.  .  FINAL. 
j,i  j.>  i> 


Maximum  integral 
max  (LANS)  =  0.19 
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Examples  of  expected  retained  fractions  across  NEDU  mesh  sizes  for  distributions  meeting 
the  overall  specifications.  The  output  file  of  this  MathCad  program  contains  46  such  distributions.  The 
Visual  Basic  program  "MeshFit"  fits  NEDU's  data  from  the  sieve  analysis  to  these  distributions.  The 
model  distribution  representing  the  best  fit  to  the  test  sample  is  then  compared  to  the  sample 
distribution.  That  comparison  is  quantified  by  assessments  of  "Goodness  of  Fit"  using  either  a  simple 
summed  error,  or  the  Chi  Square  statistical  test  (see  MeshFit  documentation). 


Definition  of  one  of  the  46  model  distributions  -  down  selected  from  1 35  candidate 
distributions. 


J  =  3 

mmMean  :=  1.3  +  (0.05-j) 

MeshMean:-  1.096+  16.9- 

(  mmMeany 

MeshMean  =  12.75 

mmMean  =  1.45 

i  :=  1 

mmSD  :=  (.05 •  i) 

mmSD  =  0.05 

8  mesh 


Fractions  retained  in  the  NEDU  sieves 
10  mesh  14  mesh 


GANS.  .  =  0 
J.i 


HANSj  .  =  0 


KANS.  .  =  1 
J.i 


30  mesh 

LANS.  .  =  1.35  x  10‘ 

j.i 


Another  model  distribution 


1 


j:=  11 

mmMean  :=  1.3  +  (0.05  j) 

mmMean  =  1 .85 

i  :=  1 

mmSD  :=  (.05-1) 

8  mesh 

10  mesh 

GANS.  . 

J.i 

=  0  HANS.  (  =  ( 

Yet  another  model  distribution 

j:=9 

mmMean  :=  1.3  +  (0.05 -j) 

mmMean  =  1.75 

i  :=  8 

mmSD  :=  (.05 -i) 

GANS.  . 
j.‘ 

=  0  HANS.  .  =  1 

j:=  6 

mmMean  :=  1.3  + (0.05 'j) 

mmMean  =  1.6 

i  :=  4 

mmSD  :=  (.05 -1) 

mmMean 


mmSD  =  0.05 

14  mesh 


KANS.  ,  =  0.99 
J.i 


i— 1 

mmMean  ) 


mmSD  =  0.4 


KANS.  .  =  0 
J.i 


^16.9 


1 


mmMean 


mmSD  =  0.2 


-  4 


GANS.  ,  =  4.83  x  10  HANS.  .  =  0.04 
J.i  j.i 


KANS.  ,  =  0.9 
J.i 


MeshMcan  =  10.23 


30  mesh 
LANSj  |  =  0 


MeshMean  =  10.75 


LANS.  .  =  0 

J.i 


MeshMean  =  11.66 


LANSj  1  =  0.07 


A-14 


j:=  11 

mmMean  :=  1.3  +  (0.05 -j) 

MeshMean  := 

1.096+  (  16.9 - ' - ] 

^  mmMean  J 

MeshMean  =  10.23 

mmMean  =  1.85 

i  :=  4 

mmSD  :=  (.05 -i) 

mmSD  =  0.2 

GANS. 

J. 

.  =  0  HANS.  .  = 

i  J.' 

0 

KANS.  .  =  0 

J.i 

LANS.  .  =  0 

j.' 

j:=9 

mmMean  :=  1.3  +  (0.05  j) 

MeshMean  := 

1.096+|  16.9 - - - ) 

^  mmMean  J 

MeshMean  =  10.75 

mmMean  =  1.75 

i  :=  4 

mmSD  :=  (.05 -i) 

mmSD  =  0.2 

GANS.  .  =  539  x  l(f 3  HANS.  .  =  0.14  KANS  =0.84  LANS  =0.01 

j,i  j.'  J.'  J’1 
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APPENDIX  B.  Normally  Distributed  4-8  Mesh  Sodalime 


Problem:  A)  Find  a  Gaussian  (normal)  distribution  that  fulfills  the  requirements  of  the  original 
specification  for  4-8  mesh  size  C02  absorbent  (sodalime). 

B)  Find  allowable  %  retention  of  absorbent  for  sieve  sizes  used  in  the  NEDU  T&E  laboratory. 

(  3.05^ 

ORIGIN®  ] 


j  :=  1 ..  15  i:=  1 ..  8  Iteration  range 

Gaussian  (normal!  distribution  parameters 
Mean  Standard  deviation 

Pj :»  3.0  +  .05- j  a  j  :=  i-.05  +  0.60 


Normal  distribution:  y(x ,  p  ,  c)  := 


yj  2-Jt-C 


(»-p) 

[2W] 


Integral  equations  relevant  to  specifications  for  10-14  mesh  absorbent 

All  mesh  sizes  expressed  in  diameter  of  mesh  screen  (in  mm). 


P  = 


3.1 
3.15 

3.2 
3.25 
3  3 
3.35 

3.4 
3.45 

3.5 
3.55 

3.6 
3.65 

3.7 
3.75. 


<0.65> 
0.7 
0.75 
0.8 
0.85 
0.9 
0.95 
,  1  , 


<4  mesh 

20 

y(*.Mj.«l)  d* 

1.76 


4-8  mesh 


8-30  mesh 


The  proposed  draft  sodalime  specification  calls  for  no  more  than  7%  larger  than  4  mesh,  no  less  than 
77%  in  the  4-8  mesh  range,  no  more  than  1 5%  between  8  and  30  mesh. 


Out  of  120  candidate  distributions,  we  found  57  that  met  all  criteria  for  the  specification,  and  were 
therefore  classified  as  model  distributions.  The  manner  in  which  these  57distributions  were  identified  are 
indicated  on  the  following  pages.  Of  those  57,  7  are  listed  on  the  last  two  pages. 
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Fraction  of  absorbent  larger  than  #4  mesh 


Matrix  of  absorbent  fraction  greater  than 
#4  mesh  (4.76  mm)  in  size  for  various  distribution 
j's  (means)  and  i's  (standard  deviations). 


4.26  x  10“  3 

7.29  x  10“  3 

0.01 

0.02 

0.02 

0.03 

0.04 

0.04 

5.33  x  10“  3 

8.86  x  10“  3 

0.01 

0.02 

0.03 

0.03 

0.04 

0.05 

6.63  x  10“  3 

0.01 

0.02 

0.02 

0.03 

0.04 

0.05 

0.05 

8.2  x  10“  3 

0.01 

0.02 

0.03 

0.03 

0.04 

0.05 

0.06 

0.01 

0.02 

0.02 

0.03 

0.04 

0.05 

0.06 

0.07 

0.01 

0.02 

0.03 

0.03 

0.04 

0.05 

0.06 

0.07 

0.02 

0.02 

0.03 

0.04 

0.05 

0.06 

0.07 

0.08 

0.02 

0.03 

0.03 

0.04 

0.05 

0.07 

0.08 

0.09 

0.02 

0.03 

0.04 

0.05 

0.06 

0.07 

0.08 

0.1 

0.03 

0.04 

0.05 

0.06 

0.07 

0.08 

0.09 

0.1 

0.03 

0.04 

0.05 

0.07 

0.08 

0.09 

0.1 

0.11 

0.04 

0.05 

0.06 

0.07 

0.09 

0.1 

0.11 

0.12 

0.04 

0.06 

0.07 

0.08 

0.1 

0.11 

0.12 

0.13 

0.05 

0.06 

0.08 

0.09 

0.11 

0.12 

0.13 

0.14 

0.06 

0.07 

0.09 

0.1 

0.12 

0.13 

0.14 

0.16 

Binary  Outcome 


f  i 

1 

1 

1 

1 

1 

1 

r 

1 

1 

l 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

0 

1 

1 

1 

1 

1 

1 

1 

0 

1 

1 

1 

1 

1 

1 

0 

0 

1 

1 

1 

1 

1 

0 

0 

0 

1 

1 

1 

1 

1 

0 

0 

0 

1 

1 

1 

1 

0 

0 

0 

0 

1 

1 

1 

0 

0 

0 

0 

0 

1 

1 

1 

0 

0 

0 

0 

0 

1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

o. 

Specific  Specification  Criterion 
XA.  .:=  if(A.  .  <  .07,  1  ,0) 


1  =  meets  specification 
0  =  does  not  meet  specification 


Each  "I"  identifies  a  theoretical  distribution  that 
meets  the  portion  of  the  specification  that 
requires  no  more  than  1%  of  the  absorbent  can 
be  larger  than  #4  mesh. 


(B1) 
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Fraction  of  absorbent  between  #4  and  #8  mesh 
(4.76  mm  to  2.36  mm) 


^0.85 

0.83 

0.81 

0.79 

0.77 

0.75 

0.73 

0.71 

0.87 

0.85 

0.82 

0.8 

0.78 

0.76 

0.74 

0.72 

0.88 

0.86 

0.84 

0.82 

0.79 

0.77 

0.75 

0.73 

0.89 

0.87 

0.85 

0.83 

0.81 

0.78 

0.76 

0.74 

0.9 

0.88 

0.86 

0.84 

0.81 

0.79 

0.77 

0.75 

0.91 

0.89 

0.87 

0.85 

0.82 

0.8 

0.78 

0.75 

0.92 

0.9 

0.88 

0.85 

0.83 

0.81 

0.78 

0.76 

0.93 

0.91 

0.88 

0.86 

0.83 

0.81 

0.79 

0.76 

0.93 

0.91 

0.89 

0.86 

0.84 

0.81 

0.79 

0.77 

0.93 

0.91 

0.89 

0.87 

0.84 

0.82 

0.79 

0.77 

0.94 

0.91 

0.89 

0.87 

0.84 

0.82 

0.79 

0.77 

0.93 

0.91 

0.89 

0.87 

0.84 

0.82 

0.79 

0.77 

0.93 

0.91 

0.89 

0.86 

0.84 

0.82 

0.79 

0.77 

0.93 

0.91 

0.88 

0.86 

0.84 

0.81 

0.79 

0.77 

^0.92 

0.9 

0.88 

0.86 

0.83 

0.81 

0.78 

0.76 

Binary  Outcome 


'] 

1 

1 

1  0 

0 

0 

<A 

1 

1 

1 

1  1 

0 

0 

0 

1 

1 

1 

1  1 

1 

0 

0 

1 

1 

1 

1  1 

1 

0 

0 

1 

1 

1 

1  1 

1 

0 

0 

1 

1 

1 

1  1 

1 

1 

0 

1 

1 

1 

1  1 

1 

1 

0 

1 

1 

1 

1  1 

1 

1 

0 

1 

1 

1 

1  1 

1 

1 

0 

1 

1 

1 

1  1 

1 

1 

0 

1 

1 

1 

1  I 

1 

1 

0 

1 

1 

1 

1  1 

1 

1 

0 

1 

1 

1 

1  1 

1 

1 

0 

1 

1 

1 

1  1 

1 

1 

0 

J 

1 

1 

1  1 

1 

1 

0, 

Specific  Specification  Criterion 
XB.  .:=  if(B.  .>0.77,1,0) 
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Fraction  of  absorbent  between  #8  and  #30  mesh 
(2.36  mm  to  0.6  mm) 


'0.14 

0.16 

0.18 

0.19 

0.21 

0.22 

0.23 

0.24 

0.13 

0.15 

0.16 

0.18 

0.19 

0.2 

0.21 

0.22 

0.11 

0.13 

0.15 

0.16 

0.17 

0.19 

0.2 

0.21 

0.1 

0.1) 

0.13 

0.15 

0.16 

0.17 

0.19 

0.2 

0.09 

0.1 

0.12 

0.13 

0.15 

0.16 

0.17 

0.18 

0.07 

0.09 

0.) 

0.12 

0.13 

0.15 

0.16 

0.17 

0.06 

0.08 

0.09 

0.11 

0.12 

0.13 

0.15 

0.16 

0.05 

0.07 

0.08 

0.1 

0.11 

0.12 

0.14 

0.15 

0.05 

0.06 

0.07 

0.09 

0.1 

0.11 

0.12 

0.14 

0.04 

0.05 

0.06 

0.08 

0.09 

0.1 

0.11 

0.13 

0.03 

0.04 

0.06 

0.07 

0.08 

0.09 

0.1 

0.12 

0.03 

0.04 

0.05 

0.06 

0.07 

0.08 

0.1 

0.11 

0.02 

0.03 

0.04 

0.05 

0.06 

0.08 

0.09 

0.1 

0.02 

0.03 

0.04 

0.05 

0.06 

0.07 

0.08 

0.09 

1,0.02 

0.02 

0.03 

0.04 

0.05 

0.06 

0.07 

0.08 

Binary  Outcome 


( 


XC  = 


0  0  0  0  0  0  o'1 
1  0  0  0  0  0  0 
1  1  0  0  0  0  0 
1  1  I  0  0  0  0 
I  1  1  I  0  0  0 
111110  0 
1  I  I  1  I  1  0 


Specific  Specification  Criterion 
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XTOT  = 


'3  2  2 
3  3  2 
3  3  3 
3  3  3 
3  3  3 
3  3  3 
3  3  3 
3  3  3 
3  3  3 
3  3  3 
3  3  3 
3  3  3 
3  3  3 
3  3  2 
^3  2  2 


2 

2 

2 

3 

3 

3 

3 

3 

3 

3 

3 

2 

2 

2 

2 


1  i  i  r 

2  111 
2  2  11 
2  2  11 
3  2  11 
3  3  2  0 
3  3  3  0 
3  3  2  1 
3  2  2  1 
3  2  2  1 
2  2  2  1 
2  2  2  1 
2  2  2  1 
2  2  2  1 

2  2  2  1  i 


FINAL.  .  :=  iff  XTOT.  .  =  3,1 
J.i  l  J.l 


.0 


I's  identify  distributions  that  simultaneously 
meet  all  three  specification  requirements. 
That  is,  these  are  good  model  distributions. 


This  matrix  is  used  below  as  a  mask  or 
filter  to  select  valid  values  for  retained 
fractions  with  NEDU's  mesh  screens. 


FINAL 


Sum  of  binary  outcomes 

XTOT.  .:=  (XA.  .  +  XB.  .  +  XC 
J.*  V  J.l  J.'  J 


/j 

0 

0 

0 

0 

0 

0 

o' 

1 

1 

0 

0 

0 

0 

0 

0 

1 

1 

1 

0 

0 

0 

0 

0 

1 

1 

1 

1 

0 

0 

0 

0 

1 

1 

1 

1 

1 

0 

0 

0 

1 

1 

1 

1 

1 

1 

0 

0 

1 

1 

1 

1 

1 

1 

1 

0 

1 

1 

1 

1 

1 

1 

0 

0 

1 

1 

1 

1 

1 

0 

0 

0 

1 

1 

1 

1 

1 

0 

0 

0 

1 

1 

I 

1 

0 

0 

0 

0 

1 

1 

1 

0 

0 

0 

0 

0 

1 

1 

l 

0 

0 

0 

0 

0 

1 

1 

0 

0 

0 

0 

0 

0 

V 

0 

0 

0 

0 

0 

0 

0, 
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The  following  definite  integrals  define  the  fraction  of  the  normal  distribution  falling  within 
the  confines  of  the  six  sieves  available  at  NEDU. 


3.5  sieve 

5  sieve 

8  sieve 

f2° 

E.  . ;=  1  y(x,Hj,<Tj)dx 
’  J5.6 

v=. 

-5.6 

y(x.Fj.Oi)  dx 

4.0 

r4.0 

Gi  i:=  y(x,Hj,crj)dx 

’  •'2.36 

10  sieve 

14  sieve 

30  sieve 

/■2.36 

H.  y(*.Fj>°i) 

’  J2.06 

V\ 

r2.06 

y(x,Hj,Oj)  dx 

1.4 

fl.4 

Lj,:=  y(*. FJ »CTi)  dx 

’  .6 

3.5  sieve  Fractions  retained  by  a  #3.5  mesh 
sieve. 


(  _  5 

4.79  x  10 

1.42  x  10"  4 

3.45  x 

io"4 

7.25  x  10"4 

1.35  x  10"  3 

2.3  x  10"  3 

3.63  x  10" 

3 

-3\ 

5.39  x  10 

6.53  x  10"  5 

1.86  x  10"  4 

4.38  x 

10"  4 

8.97  x  10“  4 

1.63  x  10"  3 

2.74  x  10"3 

4.25  x  10" 

3 

6.21  x  10"  3 

8.85  x  10"  5 

2.42  x  10"4 

5.55  x 

io- 4 

1.1  x  10"  3 

1.97  x  10"  3 

3.24  x  IO"3 

4.95  x  10" 

3 

7.14  x  10"  3 

1.19  x  10"4 

3.15  x  10"4 

6.99  x 

io- 4 

1.35  x  10"3 

2.37  x  10"  3 

3.83  x  10“  3 

5.76  x  10" 

3 

8.2  x  IO"3 

1.6  x  10"  4 

4.07  x  10"4 

8.64  x 

10"  4 

1.65  x  10"3 

2.85  x  10"  3 

4.51  x  10"  3 

6.69  x  10" 

3 

9.39  x  10"  3 

2.14  x  10“  4 

5.24  x  10"  4 

1.08  x 

io- 3 

2.02  x  10"3 

3.41  x  10"  3 

5.3  x  10"3 

7.74  x  10" 

3 

0.01 

2.83  x  10“  4 

6.71  x  10"4 

1.35  x 

io": 3 

2.46  x  10"  3 

4.06  x  10"3 

6.21  x  10"  3 

8.93  x  10" 

3 

0.01 

3.74  x  10"  4 

8.36  x  10"  4 

1.68  x 

10"  3 

2.98  x  10"  3 

4.82  x  10"  3 

7.25  x  10“  3 

0.01 

0.01 

4.91  x  10"  4 

1.07  x  10"  3 

2.07  x 

io"3 

3.6  x  10"3 

5.71  x  10"  3 

8.45  x  10"3 

0.01 

0.02 

6.41  x  10"  4 

1.35  x  10"  3 

2.55  x 

io"3 

4.33  x  10"  3 

6.74  x  10"  3 

9.82  x  10"  3 

0.01 

0.02 

8.06  x  10"  4 

1.7  x  10"3 

3.13  x 

io"3 

5.2  x  10"3 

7.94  x  10"3 

0.01 

0.02 

0.02 

1.05  x  10"3 

2.14  x  10”  3 

3.83  x 

io" 3 

6.21  x  10"  3 

9.31  x  10"  3 

0.01 

0.02 

0.02 

1.35  x  10"3 

2.67  x  10"  3 

4.66  x 

io" 3 

7.39  x  10"  3 

0.01 

0.02 

0.02 

0.03 

1.73  x  10"3 

3.32  x  10"  3 

5.65  x 

10"  3 

8.77  x  10"  3 

0.01 

0.02 

0.02 

0.03 

,2.21  x  10“  3 

4.11  x  10"  3 

6.82  x 

10" 3 

0.01 

0.01 

0.02 

0.03 

0.03  j 

B-6 


EANS  = 


Apply  binary  filter 

EANS.  .:=  E.  .-FINAL.  . 
J,i  1,1  J.' 


4.79  x  10"  5 

0 

0 

0 

0 

0 

0 

\ 

0 

6.53  x  10"  5 

1.86  x  10"  4 

0 

0 

0 

0 

0 

0 

8.85  x  10"  5 

2.42  x  10“  4 

5.55  x  10"4 

0 

0 

0 

0 

0 

1.19  x  10"  4 

3.15  x  10"4 

6.99  x  10' 4 

1.35  x  10"  3 

0 

0 

0 

0 

1.6  x  10“  4 

4.07  x  10“  4 

8.64  x  10“  4 

1.65  x  10“  3 

2.85  x  10"  3 

0 

0 

0 

2.14  x  10"  4 

5.24  x  10"  4 

1.08  x  10' 3 

2.02  x  10"  3 

3.41  x  10"  3 

5.3  x  10"3 

0 

0 

2.83  x  10”  4 

6.71  x  10"  4 

1.35  x  10“  3 

2.46  x  10"  3 

4.06  x  10"  3 

6.21  x  10"  3 

8.93  x  10"  3 

0 

3.74  x  10"4 

8.36  x  10"  4 

1.68  x  10'3 

2.98  x  10' 3 

4.82  x  10"  3 

7.25  x  10"  3 

0 

0 

4.91  x  10“  4 

1.07  x  10"3 

2.07  x  10'3 

3.6  x  10"  3 

5.71  x  10"  3 

0 

0 

0 

6.41  x  10"  4 

1.35  x  10"  3 

2.55  x  10"  3 

4.33  x  10"  3 

6.74  x  10"  3 

0 

0 

0 

8.06  x  10“  4 

1.7  x  10“  3 

3.13  x  10"3 

5.2  x  10"  3 

0 

0 

0 

0 

1.05  x  10-3 

2.14  x  10"  3 

3.83  x  10"  3 

0 

0 

0 

0 

0 

1.35  x  10"  3 

2.67  x  10"  3 

4.66  x  10"  3 

0 

0 

0 

0 

0 

1.73  x  10"  3 

3.32  x  10"  3 

0 

0 

0 

0 

0 

0 

2.21  x  10"  3 

0 

0 

0 

0 

0 

0 

0, 

B-7 


5  sieve 


'0.07 

0.09 

0.1 

0.12 

0.13 

0.14 

0.16 

0.17' 

0.08 

0.1 

0.11 

0.13 

0.14 

0.16 

0.17 

0.18 

0.1 

0.11 

0.13 

0.14 

0.16 

0.17 

0.18 

0.19 

0.11 

0.13 

0.14 

0.16 

0.17 

0.18 

0.19 

0.2 

0.12 

0.14 

0.16 

0.17 

0.19 

0.2 

0.21 

0.22 

0.14 

0.16 

0.17 

0.19 

0.2 

0.21 

0.22 

0.23 

0.16 

0.18 

0.19 

0.21 

0.22 

0.23 

0.24 

0.25 

0.18 

0.19 

0.21 

0.22 

0.24 

0.25 

0.25 

0.26 

0.2 

0.21 

0.23 

0.24 

0.25 

0.26 

0.27 

0.28 

0.22 

0.24 

0.25 

0.26 

0.27 

0.28 

0.29 

0.29 

0.24 

0.26 

0.27 

0.28 

0.29 

0.3 

0.3 

0.31 

027 

0.28 

0.29 

0.3 

0.31 

0.32 

0.32 

0.32 

0.29 

0.31 

0.32 

0.32 

0.33 

0.33 

0.34 

0.34 

0.32 

0.33 

0.34 

0.35 

0.35 

0.35 

0.35 

0.35 

,0.35 

0.36 

0.36 

0.37 

0.37 

0.37 

0.37 

0.37, 

Results  of  Binary  filtering 


'0.07 

0 

0 

0 

0 

0 

0 

o' 

0.08 

0.1 

0 

0 

0 

0 

0 

0 

Apply  binary  filter 

0.1 

0.11 

0.13 

0 

0 

0 

0 

0 

FANS.  .:=  F.  .  FINAL. 
J.i  J.i  J 

0.11 

0.13 

0.14 

0.16 

0 

0 

0 

0 

0.12 

0.14 

0.16 

0.17 

0.19 

0 

0 

0 

0.14 

0.16 

0.17 

0.19 

0.2 

0.21 

0 

0 

0.16 

0.18 

0.19 

0.21 

0.22 

0.23 

0.24 

0 

FANS  = 

0.18 

0.19 

0.21 

022 

0.24 

0.25 

0 

0 

Maximum  integral 

0.2 

0.21 

0.23 

0.24 

0.25 

0 

0 

0 

max(FANS)  =  0.35 

0.22 

0.24 

0.25 

0.26 

0.27 

0 

0 

0 

0.24 

0.26 

0.27 

0.28 

0 

0 

0 

0 

0.27 

0.28 

0.29 

0 

0 

0 

0 

0 

0.29 

0.31 

0.32 

0 

0 

0 

0 

0 

0.32 

0.33 

0 

0 

0 

0 

0 

0 

,0.35 

0 

0 

0 

0 

0 

0 

0, 
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'o^s 

0.75 

0.72 

0.69 

0.66 

0.63 

0.61 

0.58n 

0.79 

0.76 

0.72 

0.69 

0.66 

0.64 

0.61 

0.59 

0.79 

0.76 

0.73 

0.69 

0.67 

0.64 

0.61 

0.59 

0.79 

0.76 

0.73 

0.69 

0.67 

0.64 

0.61 

0.59 

0.79 

0.76 

0.72 

0.69 

0.66 

0.64 

0.61 

0.59 

0.79 

0.75 

0.72 

0.69 

0.66 

0.63 

0.61 

0.58 

0.78 

0.74 

0.71 

0.68 

0.66 

0.63 

0.6 

0.58 

0.77 

0.74 

0.71 

0.68 

0.65 

0.62 

0.6 

0.58 

0.75 

0.72 

0.7 

0.67 

0.64 

0.62 

0.59 

0.57 

0.74 

0.71 

0.68 

0.66 

0.63 

0.61 

0.59 

0.56 

0.72 

0.7 

0.67 

0.64 

0.62 

0.6 

0.58 

0.56 

0.7 

0.68 

0.65 

0.63 

0.61 

0.59 

0.57 

0.55 

0.68 

0.66 

0.64 

0.62 

0.6 

0.58 

0.56 

0.54 

0.66 

0.64 

0.62 

0.6 

0.58 

0.56 

0.54 

0.53 

v0.63 

0.62 

0.6 

0.58 

0.56 

0.55 

0.53 

0.52, 

0.78 

0 

0 

0 

0 

0 

0 

0^ 

0.79 

0.76 

0 

0 

0 

0 

0 

0 

0.79 

0.76 

0.73 

0 

0 

0 

0 

0 

0.79 

0.76 

0.73 

0.69 

0 

0 

0 

0 

0.79 

0.76 

0.72 

0.69 

0.66 

0 

0 

0 

0.79 

0.75 

0.72 

0.69 

0.66 

0.63 

0 

0 

0.78 

0.74 

0.71 

0.68 

0.66 

0.63 

0.6 

0 

9S  = 

0.77 

0.74 

0.71 

0.68 

0.65 

0.62 

0 

0 

0.75 

0.72 

0.7 

0.67 

0.64 

0 

0 

0 

0.74 

0.71 

0.68 

0.66 

0.63 

0 

0 

0 

0.72 

0.7 

0.67 

0.64 

0 

0 

0 

0 

0.7 

0.68 

0.65 

0 

0 

0 

0 

0 

0.68 

0.66 

0.64 

0 

0 

0 

0 

0 

0.66 

0.64 

0 

0 

0 

0 

0 

0 

^0.63 

0 

0 

0 

0 

0 

0 

0, 

Apply  binary  filter 

GANS.  G.  .  FINAL. 
J,i  J.i  J. 


Maximum  integral 
max  (GANS)  =0.79 
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10  sieve 


'0.08 

0.08 

0.09 

0.09 

0.09 

0.09 

0.09 

0.08 ' 

0.07 

0.08 

0.08 

0.08 

0.08 

0.08 

0.08 

0.08 

0.07 

0.07 

0.07 

0.08 

0.08 

0.08 

0.08 

0.08 

0.06 

0.06 

0.07 

0.07 

0.07 

0.07 

0.07 

0.07 

0.05 

0.06 

0.06 

0.06 

0.07 

0.07 

0.07 

0.07 

0.05 

0.05 

0.06 

0.06 

0.06 

0.06 

0.07 

0.07 

0.04 

0.05 

0.05 

0.05 

0.06 

0.06 

0.06 

0.06 

0.04 

0.04 

0.05 

0.05 

0.05 

0.06 

0.06 

0.06 

0.03 

0.04 

0.04 

0.05 

0.05 

0.05 

0.05 

0.06 

0.03 

0.03 

0.04 

0.04 

0.04 

0.05 

0.05 

0.05 

0.02 

0.03 

0.03 

0.04 

0.04 

0.04 

0.05 

0.05 

0.02 

0.02 

0.03 

0.03 

0.04 

0.04 

0.04 

0.05 

0.02 

0.02 

0.03 

0.03 

0.03 

0.04 

0.04 

0.04 

0.01 

0.02 

0.02 

0.03 

0.03 

0.03 

0.04 

0.04 

<0.01 

0.02 

0.02 

0.02 

0.03 

0.03 

0.03 

0.04, 

HANS  = 


'0.08 

0 

0 

0 

0 

0 

0 

o' 

0.07 

0.08 

0 

0 

0 

0 

0 

0 

0.07 

0.07 

0.07 

0 

0 

0 

0 

0 

0.06 

0.06 

0.07 

0.07 

0 

0 

0 

0 

0.05 

0.06 

0.06 

0.06 

0.07 

0 

0 

0 

0.05 

0.05 

0.06 

0.06 

0.06 

0.06 

0 

0 

0.04 

0.05 

0.05 

0.05 

0.06 

0.06 

0.06 

0 

0.04 

0.04 

0.05 

0.05 

0.05 

0.06 

0 

0 

0.03 

0.04 

0.04 

0.05 

0.05 

0 

0 

0 

0.03 

0.03 

0.04 

0.04 

0.04 

0 

0 

0 

0.02 

0.03 

0.03 

0.04 

0 

0 

0 

0 

0.02 

0.02 

0.03 

0 

0 

0 

0 

0 

0.02 

0.02 

0.03 

0 

0 

0 

0 

0 

0.01 

0.02 

0 

0 

0 

0 

0 

0 

,0.01 

0 

0 

0 

0 

0 

0 

o. 

Apply  binary  filter 

HANS.  .:=  H.  .  FINAL.  . 
j.l  j.<  j.i 


Maximum  integral 

max(HANS)  =  0.08 


B-10 


14  sieve 


KANS  = 


f  0.06 

0.07 

0.08 

0.09 

0.1 

0.1 

0.11 

0.11 ' 

0.05 

0.06 

0.07 

0.08 

0.09 

0.09 

0.1 

0.1 

0.04 

0.05 

0.06 

0.07 

0.08 

0.09 

0.09 

0.1 

0.04 

0.05 

0.06 

0.06 

0.07 

0.08 

0.09 

0.09 

0.03 

0.04 

0.05 

0.06 

0.07 

0.07 

0.08 

0.08 

0.03 

0.03 

0.04 

0.05 

0.06 

0.07 

0.07 

0.08 

0.02 

0.03 

0.04 

0.05 

0.05 

0.06 

0.07 

0.07 

0.02 

0.03 

0.03 

0.04 

0.05 

0.06 

0.06 

0.07 

K  = 

0.02 

0.02 

0.03 

0.04 

0.04 

0.05 

0.06 

0.06 

0.01 

0.02 

0.02 

0.03 

0.04 

0.04 

0.05 

0.06 

0.01 

0.02 

0.02 

0.03 

0.03 

0.04 

0.05 

0.05 

8.56  x  10“ 

3 

0.01 

0.02 

0.02 

0.03 

0.04 

0.04 

0.05 

6.95  x  10" 

3 

0.01 

0.02 

0.02 

0.03 

0.03 

0.04 

0.04 

3 

-  3 

5.62  x  10“ 

9.06  x  10 

0.01 

0.02 

0.02 

0.03 

0.03 

0.04 

^4.51  x  10' 

3 

7.49  x  10“  3 

0.01 

0.02 

0.02 

0.03 

0.03 

0.04  j 

0.06 

0 

0 

0 

0 

0 

0 

0> 

0.05 

0.06 

0 

0 

0 

0 

0 

0 

0.04 

0.05 

0.06 

0 

0 

0 

0 

0 

0.04 

0.05 

0.06 

0.06 

0 

0 

0 

0 

Apply  binary  filter 

0.03 

0.04 

0.05 

0.06 

0.07 

0 

0 

0 

KANS.  .  :=  K.  .  FINAL. 

J.i  J.i  J 

0.03 

0.03 

0.04 

0.05 

0.06 

0.07 

0 

0 

0.02 

0.03 

0.04 

0.05 

0.05 

0.06 

0.07 

0 

0.02 

0.03 

0.03 

0.04 

0.05 

0.06 

0 

0 

0.02 

0.02 

0.03 

0.04 

0.04 

0 

0 

0 

0.01 

0.02 

0.02 

0.03 

0.04 

0 

0 

0 

0.01 

0.02 

0.02 

0.03 

0 

0 

0 

0 

Maximum  integral 

8.56  x  10“  3 

0.01 

0.02 

0 

0 

0 

0 

0 

raax(KANS)  =  0.07 

6.95  x  10“  3 

0.01 

0.02 

0 

0 

0 

0 

0 

_  3 

-3 

5.62  x  10 

9.06  x  10 

0 

0 

0 

0 

0 

0 

_  3 

4.51  x  10 

0 

0 

0 

0 

0 

0 

0. 

B-1 1 


30  sieve 


(  -3 

5.49  x  10 

8.98  x  10"  3 

0.01 

0.02 

0.02 

0.03 

0.04 

\ 

0.04 

4.4  x  10"  3 

7.4  x  10"  3 

0.01 

0.02 

0.02 

0.03 

0.03 

0.04 

3.5  x  10"  3 

6.07  x  10“  3 

9.48  x  10"  3 

0.01 

0.02 

0.02 

0.03 

0.03 

2.78  x  10”3 

4.96  x  10"3 

7.93  x  10"  3 

0.01 

0.02 

0.02 

0.03 

0.03 

2.19  x  10"  3 

4.03  x  10"  3 

6.61  x  10"  3 

9.91  x  10"3 

0.01 

0.02 

0.02 

0.03 

1.72  x  10"3 

3.26  x  10"  3 

5.49  x  10"  3 

8.41  x  10"  3 

0.01 

0.02 

0.02 

0.03 

1.34  x  10"  3 

2.63  x  10"  3 

4.54  x  10"3 

7.1  x  10"  3 

0.01 

0.01 

0.02 

0.02 

L  = 

1.04  x  10"3 

2.11  x  10“  3 

3.74  x  I0"3 

5.98  x  10"  3 

8.82  x  10”  3 

0.01 

0.02 

0.02 

8  x  10" 4 

1.68  x  10"  3 

3.06  x  10"  3 

5.01  x  10"  3 

7.54  x  10"  3 

0.01 

0.01 

0.02 

6.13  x  10"4 

1.33  x  10"3 

2.5  x  10"  3 

4  .19  x  10”  3 

6.42  x  10"  3 

9  .18  x  10"  3 

0.01 

0.02 

4.68  x  10"4 

1.05  x  10"  3 

2.03  x  10"  3 

3.49  x  10"3 

5.45  x  10"  3 

7.93  x  10”  3 

0.01 

0.01 

3.54  x  10"  4 

8.27  x  10“  4 

1.65  x  10”  3 

2.89  x  10"3 

4.62  x  10"  3 

6.82  x  10"3 

9.49  x  10" 

3 

0.01 

2.67  x  10"  4 

6.47  x  10"  4 

1.33  x  10”  3 

2.39  x  10"  3 

3.89  x  10"  3 

5.86  x  10”  3 

8.27  x  10” 

3 

0.01 

-  4 

2x  10 

5.04  x  10"4 

1.06  x  10"  3 

1.97  x  10"3 

3.27  x  10”  3 

5.01  x  10"  3 

7.19  x  10" 

3 

9.76  x  10"  3 

^  1 .49  x  10"4 

3.9  x  10"  4 

8.51  x  10“  4 

1.61  x  10"3 

2.74  x  10"  3 

4.28  x  10”  3 

6.23  x  10" 

3 

8.57  x  10”  3 , 

B-12 


LANS 


Apply  binary  filter 

LANS.  .:=  L.  FINAL.  . 
j,i  j,i  J.> 

(  -  3 

5.49  x  10  0  0 


4.4  x  10  3 

7.4  x 

10  3 

0 

3.5  x  10"  3 

6.07  x 

io- : 3 

9.48  x  10"  3 

2.78  x  10”  3 

4.96  x 

id"  ; 3 

7.93  x  IO"3 

2.19  x  10"  3 

X 

© 

rr 

10"  3 

6.61  x  10"  3 

1.72  x  10-3 

3.26  x 

10"  3 

5.49  x  10"  3 

1.34  x  10"3 

2.63  x 

10"  3 

4.54  x  IO"3 

1.04  x  10"3 

2.11  x 

10"  3 

3.74  x  IO"3 

8  x  10“ 4 

1.68  x 

10"  3 

3.06  x  10"  3 

6.13  x  10"4 

1.33  x 

,o"3 

2.5 x  IO"3 

4.68  x  10"  4 

© 

yj\ 

X 

io" 3 

2.03  x  IO"3 

3.54  x  10"  4 

8.27  x 

io"4 

1.65  x  10"  3 

2.67  x  10"4 

6.47  x 

10- 4 

1.33  x  10"  3 

2x  10"  4 

5.04  x 

io" 4 

0 

J.49x  10"4 

0 

0 

Maximum  integral 
max(LANS)  =0.02 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0.0 1 

0 

0 

0 

0 

9.91  x  10"  3 

0.01 

0 

0 

0 

8.41  x  10"  3 

0.01 

0.02 

0 

0 

7.1  x  10"  3 

0.01 

0.01 

0.02 

0 

5.98  x  10"  3 

8.82  x  10" 

3 

0.01 

0 

0 

5.01  x  10"  3 

7.54  x  10" 

3 

0 

0 

0 

4.19  x  10"  3 

6.42  x  10" 

3 

0 

0 

0 

3.49  x  10"  3 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0> 
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Examples  of  expected  retained  fractions  across  NEDU  mesh  sizes  for  distributions  meeting 
the  overall  specifications.  The  output  file  of  this  MathCad  program  contains  57  such  distributions.  The 
Visual  Basic  program  "MeshFit"  fits  NEDU's  data  from  the  sieve  analysis  to  these  distributions.  The 
model  distribution  representing  the  best  fit  to  the  test  sample  is  then  compared  to  the  sample 
distribution.  That  comparison  is  quantified  by  assessments  of  "Goodness  of  Fit"  using  either  a  simple 
summed  error,  or  the  Chi  Square  statistical  test  (see  MeshFit  documentation). 


Definition  of  seven  of  the  57  model  distributions  -  down  selected  from  120  candidate 
distributions. 


Example  B1 


j  :=  1  mmMean  :=  2.95  +  (0.065  j)  MeshMean  :=  1.096  +  ^16.9- 

mmMean  =  3.02 

i  :=  1  mmSD  :=  0.05 -i  +  0.60  mmSD  =  0.65 

Fractions  retained  in  the  NEDU  sieves 

3.5  mesh 

5  mesh 

8  mesh 

EANS.  .  =  4.79  x  10~5 

FANS.  .  =  0.07 

GANS.  . 

=  0.78 

J.i 

J.i 

J.' 

Example  B2 

j  :=  7  mmMean  :=  2.95  +  (0.065-j) 

MeshMean  := 

1.096  + 

mmMean  = 

=  3.41 

i  :=  1  mmSD  :=  0.05  1  +  0.60 

mmSD  =  0.65 

3.5  mesh 

5  mesh 

8  mesh 

EANS.  .  =  2.83  x  10”  4 

FANS.  .  =  0.16 

GANS.  . 

=  0.78 

J.i 

J.i 

J.i 

Example  B3 

j  :=  15  mmMean 

:=  2.95  +  (0.065  j) 

MeshMean  := 

1 .096  ■ 

mmMean 

=  3.93 

i  :=  1  mmSD  := 

0.05 -i  +  0.60 

mmSD  =  0.65 

3.5  mesh 

5  mesh 

8  mesh 

EANS.  .  =  2.21  x  10' 3 

FANS.  .  =  0.35 

GANS.  ,= 

0.63 

J.i 

J.' 

J.i 

I 


mmMean  j 


MeshMean  =  6.7 


1 0  mesh 

HANS.  .  =  0.08 

J.  • 


1 4  mesh 

KANS.  .  =  0.06 

J.i 


1 


mmMean , 


MeshMean  =  6.06 


10  mesh 
HANS.  .  =  0.04 


14  mesh 
KANS.  .  =  0.02 


f  16-9 - ! - 1 

V,  mmMean  J 


1 0  mesh 


MeshMean  =  5.4 


14  mesh 


HANS.  .  =  0.01  KANS.  .  =  4.51  x  10 

l.i  J.i 


-3 
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Example  B4 


j  :=  4  mmMean  :=  2.95  +  (0.065-j) 

mmMean  =  3.21 

i  :=  4  ramSD  :=  0.05-1  +  0.60 


3.5  mesh 

EANS.  .=  1.35  x  10" 

J.' 


5  mesh 

FANS.  .  =  0.16 

J.' 


MeshMean  :=  1 .096  + 


1 


mmMean 


MeshMean  =  6.36 


mmSD  =  0.8 
8  mesh 

GANS.  .  =  0.69 
J.i 


1 0  mesh 


14  mesh 


HANS.  .  =  0.07 
J.i 


KANS.  .  =  0.06 


Example  B5 

j  :=  7  mmMean  :=  2.95  +  (0.065-j) 

mmMean  =  3.41 

i  :=  4  mmSD  :=  0.05 -i  +  0.60 


MeshMean  :=  1.096  + 


1 


mmMean 


MeshMean  =  6.06 


mmSD  =  0.8 


3.5  mesh 

5  mesh 

EANS.  .  =  2.46  x  10-3 

FANS.  .  =  0.21 

8  mesh 


10  mesh 


14  mesh 


GANS.  ,  =  0.68 
J.i 


HANS.  .  =  0.05 

J.' 


KANS.  j  =  0.05 


Example  B6 


j:=  11 

mmMean  :=  2.95  +  (0.065-j) 

MeshMean  :=  1.096  +  | 

f 

16.9- 

1  ") 

MeshMean  =  5.71 

mmMean  ) 

mmMean  =  3.67 

< 

i  :=  4 

mmSD:=  0.05-1  +  0.60 

mmSD  =  0.8 

3.5  mesh 

5  mesh 

8  mesh 

1 0  mesh 

14  mesh 

EANS.  .  =  5.2  x  10"  3  FANS.  .  =  0.28 

J.i  J.' 

GANS.  ,  =  0.64 

J.' 

HANS.  .  =  0.04 
J.' 

KANS.  .  =  0.03 

J.i 

Example  B7 

mmMean  :=  2.95  +  (0.065-j) 

mmMean  =  3.41 

MeshMean  :=  1.096  +  j 

f 

16.9 

< 

1  } 

MeshMean  =  6.06 

j  :=  7 

mmMean  ) 

i  :=  7 

mmSD  :=  0.05  -i  +  0.60 

mmSD  =  0.95 

3.5  mesh 

5  mesh 

8  mesh 

10  mesh 

14  mesh 

EANS.  .  =  8.93 

x  10-3  FANS.  .  =  0.24 

J.i 

GANS.  .  =  0.6 

J.i 

HANS.  .  =  0.06 
J.' 

KANS.  .  =  0.07 

J.' 
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APPENDIX  C.  LogNormally  Distributed  10-14  Mesh  Sodalime 

Problem:  A)  Find  a  Log-Gaussian  (log-normal)  distribution  that  fulfills  the  requirements  of  the  original 
specification  for  8-12  (10-14  US)  mesh  size  C02  absorbent  (sodalime). 

B)  Find  allowable  %  retention  of  absorbent  for  sieve  sizes  used  in  the  NEDU  T&E  laboratory. 


ORIGINS  1 


j  :=  I ..  12  i  :=  1 ..  8  Iteration  range 

Log-Gaussian  (log-normal)  distribution  parameters 
Mean  Standard  deviation 

|ij  :=  .32  +  .03-j  Oj:=  1  0.025-.02 


Log-Normal  distribution: 


(ln(*)-n)2 


Integral  equations  relevant  to  specifications  for  10-14  mesh  absorbent 

All  mesh  sizes  expressed  in  diameter  of  mesh  screen  (in  mm). 


'0.35' 

-  3  ^ 

0.38 

5  x  10 

0.41 

0.03 

0.44 

0.06 

0.47 

0.08 

0.5 

0.11 

n  = 

0.53 

0.13 

0.56 

0.16 

0.59 

0.18  j 

0.62 

0.65 

,0.68, 

<7  mesh 

7-10  mesh 

10-14  mesh 

(-20 

(•28 

f2.0 

Ai  i:*  y(*.Rj.°i)d* 

Bi  i:s  y(*»Pj.°i) d* 

Ci,H 

J2.80 

J2.0 

y1.4 

14-40  mesh 
fl.4 

D:  i:=  y(x,nj  .CTj)  dx 

J  0.425 


The  proposed  draft  sodalime  specification  calls  for  no  more  than  1%  larger  than  7  mesh,  no  more  than 
30%  in  the  7-10  mesh  size  range,  no  less  than  48%  in  the  10-14  range,  no  more  than  20%  between  14 
and  40  mesh. 

Out  of  135  candidate  distributions,  we  found  46  that  met  all  criteria  for  the  specification,  and  were 
therefore  classified  as  model  distributions  The  manner  in  which  these  46  distributions  were  identified  are 
indicated  on  the  following  pages.  Of  those  46,  6  are  listed  on  the  last  two  pages. 
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Fraction  of  absorbent  larger  than  #7  mesh 


Matrix  of  absorbent  fraction  greater  than 
#7  mesh  (2.8  mm)  in  size  for  various  distribution 
j's  (means)  and  i's  (standard  deviations). 


( 

0 

0 

0 

0 

4.97  x  10"  10 

5.85  x 

.o'7 

2.87  x  10"  5 

-4\ 

3.13  x  10 

0 

0 

0 

3.92  x 

lO"15 

3.04  x  10"  9 

1.9  x 

io-6 

6.61  x  10"  5 

4.6  x  10"4 

0 

0 

0 

7.68  x 

io- 14 

1.71  x  10"8 

5.88  x 

10- 6 

1.47  x  10"4 

8.61  x  10"  4 

0 

0 

0 

-  12 

1.31  x  10 

8.84  x  IO-8 

1.72  x 

io"5 

3.14  x  IO-4 

1.55  x  10"  3 

0 

0 

0 

1.93  x 

io-" 

4.22  x  10"  7 

4.78  x 

10- 5 

4.64  x  IO-4 

2.77  x  IO-3 

0 

0 

0 

2.48  x 

10" 10 

1.86  x  IO-6 

1.26  x 

io" 4 

9.24  x  10"  4 

4.81  x  IO-3 

A  — 

0 

0 

0 

2.77  x 

io" 9 

7.54  x  IO-6 

3.15  x 

10"  4 

1.85  x  10"3 

8.15  x  10"3 

0 

0 

0 

2.69  x 

10" 8 

2.82  x  10"5 

4.78  x 

10- 4 

3.58  x  10"3 

0.01 

0 

0 

-  14 

1.59  X  10 

2.26  x 

10"  7 

9.72  x  10"  5 

1.05  x 

10- 3 

6.7  x  10"  3 

0.02 

0 

0 

1.07  x  10"  12 

1.66  x 

10- 6 

3.09  x  10-4 

2.36  x 

10- 3 

0.01 

0.03 

0 

0 

5.37  x  10"  11 

1.05  x 

io" 5 

5.14  x  10"  4 

5.09  x 

10- 3 

0.02 

0.05 

<0 

0 

-  o 

2  x  10 

5.82  x 

10- 5 

1.25  x  10"  3 

0.0 1 

0.04 

0.08  , 

Binary  Outcome 


Specific  Specification  Criterion 


xa.  . : 

J.i 


/| 

1 

1 

1 

1 

1 

1 

1  \ 

1 

1 

1 

1 

1 

1 

1 

1 

I 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

0 

1 

1 

1 

1 

1 

1 

1 

0 

1 

1 

1 

1 

1 

1 

0 

0 

1 

1 

1 

1 

1 

1 

0 

0 

V  ‘ 

1 

1 

1 

1 

0 

0 

0; 

1  =  meets  specification 
0  =  does  not  meet  specification 


Each  "1"  identifies  a  theoretical  distribution  that 
meets  the  portion  of  the  specification  that 
requires  no  more  than  1%  of  the  absorbent  can 
be  larger  than  #7  mesh. 


C-2 


Fraction  of  absorbent  between  #7  and  #10  mesh 
(2.8  mm  to  2.0  mm) 


/ 

0 

0 

4.46  x  10“  10 

1.82  x  10“ 

5 

1.11  x  10“ 

3 

8.65  x  10“  3 

0.03 

\ 

0.06 

0 

0 

1.26  x  10“  8 

9.23  x  10“ 

5 

2.95  x  10“ 

3 

0.02 

0.05 

0.09 

0 

0 

2.66  x  10“  7 

4.09  x  10“ 

4 

7.24  x  10" 

3 

0.03 

0.07 

0.12 

0 

0 

4.22  x  10“  6 

1.59  x  10“ 

3 

0.02 

0.05 

0.11 

0.17 

0 

1.18  x  10“  13 

5.03  x  10“  5 

5.41  x  10“ 

3 

0.03 

0.09 

0.16 

0.23 

B  = 

0 

1.32  x  10“  10 

4.51  x  10“  4 

0.02 

0.07 

0.15 

0.23 

031 

0 

5.66  x  10“  8 

3.06  x  10“  3 

0.04 

0.13 

0.22 

0.32 

0.4 

0 

9.28  x  10“  6 

0.02 

0.1 

0.22 

0.33 

0.42 

0.5 

0 

5.92  x  10“  4 

0.06 

0.2 

0.34 

0.46 

0.55 

0.61 

0 

0.01 

0.19 

0.38 

0.52 

0.62 

0.69 

0.74 

0 

0.15 

0.45 

0.62 

0.73 

0.81 

0.85 

0.87 

<8.57 

x  10“  3 

0.67 

0.84 

0.92 

0.98 

1.01 

1.02 

1.01 , 

Binary  Outcome 


/j 

1 

1 

1 

1 

1 

1 

r 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

0 

1 

1 

1 

1 

1 

1 

0 

0 

1 

1 

1 

1 

1 

0 

0 

0 

1 

1 

1 

1 

0 

0 

0 

0 

1 

1 

1 

0 

0 

0 

0 

0 

1 

1 

0 

0 

0 

0 

0 

0 

\1 

0 

0 

0 

0 

0 

0 

o, 

Specific  Specification  Criterion 


C-3 


Fraction  of  absorbent  between  #10  and  #14  mesh 
(2.0  mm  to  1.4  mm) 


1.41 

0.97 

0.88 

0.85 

0.84 

0.84 

0.83 

0.81  N 

1.46 

1.36 

1.17 

1.08 

1.02 

0.98 

0.94 

0.9 

1.51 

1.5 

1.39 

1.27 

1.19 

1.12 

1.05 

0.98 

1.55 

1.55 

1.51 

1.42 

1.33 

1.23 

1.14 

1.05 

1.6 

1.6 

1.59 

1.54 

1.44 

1.32 

1.21 

1.1 

1.65 

1.65 

1.65 

1.61 

1.51 

1.38 

1.25 

1.13 

1.7 

1.7 

1.7 

1.65 

1.54 

1.4 

1.26 

1.14 

1.75 

1.75 

1.74 

1.65 

1.52 

1.38 

1.25 

1.13 

1.8 

1.8 

1.74 

1.6 

1.46 

1.32 

1.2 

1.09 

1.86 

1.84 

1.67 

1.49 

1.35 

1.23 

1.13 

1.04 

1.92 

1.76 

1.47 

1.3 

1.19 

1.11 

1.04 

0.97 

1.97 

1.3 

1.13 

1.06 

1.01 

0.97 

0.93 

0.88  , 

Binary  Outcome 


/i 

1 

1 

1 

1 

1 

1 

.  \ 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

I 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

I 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

J 

1 

1 

1 

1 

1 

1 

1  / 

Specific  Specification  Criterion 


C-4 


Fraction  of  absorbent  between  #14  and  #30  mesh 
(1.40  mm  to  0.60  mm) 


(  A  11  ..  ^  /.  ..  A  CA  A  Cl  A  CO  A  CO  A  CO  A  CO  1 


4.77  x  10 

0.45 

0.54 

0.57 

0.58 

0.58 

0.58 

0.58 

0 

0.1 

0.29 

0.39 

0.44 

0.47 

0.49 

0.5 

0 

9.88  x  10“  3 

0.12 

0.24 

0.32 

0.37 

0.4 

0.43 

0 

4.14  x  10“  4 

0.04 

0.13 

0.21 

0.28 

0.32 

0.36 

0 

7.02  x  10“  6 

0.01 

0.06 

0.14 

0.2 

0.25 

0.29 

0 

4.65  x  10“  8 

2.03  x  10“  3 

0.03 

0.08 

0.14 

0.19 

0.23 

D  = 

0 

1.17  x  10"  10 

3  x  10“  4 

0.01 

0.04 

0.09 

0.14 

0.18 

0 

1.1  x  10“  13 

3.36  x  10“  5 

3.56  x  10”3 

0.02 

0.06 

0.1 

0.14 

0 

0 

2.84  x  10”  6 

1.05  x  10“  3 

0.01 

0.03 

0.07 

0.1 

0 

0 

1.81  x  10“  7 

2.7  x  10“  4 

4.7  x  10”  3 

0.02 

0.04 

0.07 

0 

0 

8.65  x  10“  9 

6.11  x  10“  5 

1.92  x  10”  3 

0.01 

0.03 

0.05 

^  o 

0 

3.12  x  10”  10 

1.21  x  10“  5 

7.28  x  10“  4 

5.54  x  10“  3 

0.02 

0.04 

Binary  Outcome 


( | 

0 

0 

0 

0 

0 

0 

o' 

1 

1 

0 

0 

0 

0 

0 

0 

1 

1 

1 

0 

0 

0 

0 

0 

1 

1 

I 

1 

0 

0 

0 

0 

1 

1 

1 

1 

1 

1 

0 

0 

1 

1 

1 

1 

1 

1 

1 

0 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

I 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1  y 

Specific  Specification  Criterion 

xd.  .:=  lf(D.  .<  .20, 1 ,0) 

J.'  V  J.»  I 


C-5 


Sum  of  binary  outcomes 


xtot  = 


r4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 


14 


3  3  3  3  3 

4  3  3  3  3 
4  4  3  3  3 
4  4  4  3  3 
4  4  4  4  4 
4  4  4  4  4 
4  4  4  4  4 
4  4  4  4  3 
4  4  4  3  3 
4  4  3  3  3 
4  3  3  3  3 
3  3  3  3  2 


3  3 ' 
3  3 
3  3 
3  3 

3  3 

4  2 
3  3 
3  2 
3  2 
2  2 
2  2 
2  2, 


xtot.  . :=  /xa.  .  +  xb.  .  +  xc.  . 
J.i  l  J.i  J.i  J.I 


Final  .  :=  if^xtot.  .  c  4 , 1  ,  Oj 


1's  identify  distributions  that  simultaneously 
meet  all  four  specification  requirements. 
That  is,  these  are  good  model  distributions. 

This  matrix  is  used  below  as  a  mask  or 
filter  to  select  valid  values  for  retained 
fractions  with  NEDU’s  mesh  screens. 


Final  - 


/. 

0 

0 

0 

0 

0 

0 

0^ 

1 

l 

0 

0 

0 

0 

0 

0 

1 

l 

1 

0 

0 

0 

0 

0 

1 

t 

1 

1 

0 

0 

0 

0 

1 

i 

1 

1 

] 

1 

0 

0 

1 

l 

1 

1 

1 

1 

1 

0 

1 

l 

1 

1 

1 

1 

0 

0 

1 

l 

1 

1 

1 

0 

0 

0 

1 

l 

1 

1 

0 

0 

0 

0 

1 

l 

1 

0 

0 

0 

0 

0 

1 

l 

0 

0 

0 

0 

0 

0 

J 

0 

0 

0 

0 

0 

0 

o. 

+  xd 


C-6 


The  following  definite  integrals  define  the  fraction  of  the  normal  distribution  falling  within 
the  confines  of  the  six  sieves  available  at  NEDU. 


3.5  sieve 

5  sieve 

8  sieve 

(•20 

r5.6 

r4.0 

y(x,nj,Oi)  dx 

5.6 

F:  y(*»f>j.°i)d* 

’  J4.0 

Gj  i:=  y(x,jij,Oi)dx 

J2.36 

10  sieve 

14  sieve 

30  sieve 

r2.36 

r206 

fl.4 

H. 

y(x,Hj,Oi)  dx 

2.06 

K :  y(x,Hj,Oj)dx 

’  J1.4 

Lj  .:=J  y(x,Mj»crj)  dx 
.6 

C-7 


3.5  sieve 


0  0  0  0  0 

0  0  0  0  0 

0  0  0  0  0 

0  0  0  0  0 

0  0  0  0  0 

0  0  0  0  0 

0  0  0  0  0 

0  0  0  0  0 

0  0  0  0  0 

0  0  0  0  0 

0  0  0  0  0 

,0  0  0  0  0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

4.65  X  10"  15 


0 


0 

0 


0 

2.48  x  10' 
1.17  x  10' 


5.33  x  10 

2.34  x  10' 


9.9  x  10 
4.04  x  10' 


1.59  x  10 
6.02  x  10 


8.31  x  10 


14 


2.95  x  10 


13 


1.02  x  10 


12 


3.44  x  10 


-  12 


1.13  x  10 


11 


3.59  x  10 


11 


1.12  x  10 


10 


3.37  x  10 


-  10 


9.92  x  )0 


10 


2.84  x  10 


-9 


7.92  x  10"  9 
2.15  x  10“  % 


ractions  retained  by  a  #3.5  mesh 
ieve. 


Results  of  binary  filtering 


Apply  binary  filter 

cans.  .  :=  E.  .  Final.  . 
J.i  J.i  J.l 


'0  0  0  0  0  0 
0  0  0  0  0  0 
0  0  0  0  0  0 
0  0  0  0  0  0 
0  0  0  0  0  0 


0  o' 
0  0 
0  0 
0  0 
0  0 


0  0  0 

0  0  0 

0  0  0 

0  0  0 

0  0  0 

,0  0  0 


0 

0 

0 

0 

0 

0 

0 


0  0  1.17  X  10"  14  0 
0  0  0  0 

0  0  0  0 

0  0  0  0 

0  0  0  0 

0  0  0  0 

0  0  0  0, 


Write  results  to  a  data  file  for 
use  by  MeshFit  software. 

\VRITEPRN("meLnl014.pm")  :=  cans 


C-8 


5  sieve 


( 

0 

0 

0 

0 

0 

3.19  x  10  15 

4.7  x 

o"11 

_  g  \ 

1.76  x  10 

0 

0 

0 

0 

0 

-  ]4 

2.01  x  10 

1.73  x 

10"  10 

4.67  x  10"  8 

0 

0 

0 

0 

0 

1.2  x  10"  13 

6.15  x 

io"10 

1.2 x  IO"7 

0 

0 

0 

0 

0 

6.83  x  10"  13 

2.11  x 

io"9 

3.02  x  10"  7 

0 

0 

0 

0 

0 

3.68  x  10"  12 

6.95  x  10  9 

7.39  x  10"  7 

F  = 

0 

0 

0 

0 

0 

1.89  x  10"  11 

2.21  x 

io"8 

1.76  x  10"  6 

0 

0 

0 

0 

0 

9.16  x  10"  11 

6.79  x 

10"  8 

4.07  x  10"  6 

0 

0 

0 

0 

722  x  10" 

15 

4.22  x  10"  10 

2.01  x 

io"7 

9.17  x  10"  6 

0 

0 

0 

0 

6.81 

x  10" 

14 

1.85  x  10"  9 

5.73  x 

io"7 

2.01  x  10“  5 

0 

0 

0 

0 

5.92 

x  10" 

13 

7.67  x  10"  9 

1.58  x 

io"6 

4.31  x  10”  5 

0 

0 

0 

0 

4.76  x  10" 

12 

3.03  x  10"  8 

4.19  x 

io" 6 

8.96  x  10"  5 

<0 

0 

0 

0 

3.53 

x  10 

11 

1.13  x  10"  7 

1.07  x 

10"  5 

1.82  x  10"  4  j 

Results  of  Binary  filtering 

Apply  binary  filter 

r0 

0 

0 

0 

0 

0 

0 

O'1 

fans.  .  :=  F.  .-Final.  . 

J.i  J.i  j.i 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

3.68  x  10"  12 

0 

0 

Maximum  integral 

0 

0 

0 

0 

0 

1.89  x  10“  11 

2.21  x  IO"8 

0 

max(fans)  =  2.21  x 

0 

0 

0 

0 

0 

9.16  x  10"M 

0 

0 

0 

0 

0 

0 

7.22  x  10"  15 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

<0 

0 

0 

0 

0 

0 

0 

0, 

C-9 


8  sieve 


( 

0 

0 

0 

2.47  x  10”  10 

1.53  x  I0”6 

1.11  x  10”  4 

1.27  x  10”  3 

-3Nt 

5.88  x  10 

0 

0 

0 

2.64  x  10”9 

6.21  x  10”6 

2.82  x  10”4 

2.48  x  10”  3 

9.78  x  10”  3 

0 

0 

0 

2.46  x  10“  8 

233  x  10”  5 

6.81  x  10”  4 

4.69  x  10”  3 

0.02 

0 

0 

-  14 

3.64  x  10 

2  x  10"7 

8.08  x  10”  5 

1.56  x  10”  3 

8.56  x  10”  3 

0.03 

0 

0 

2.09  x  10”  12 

1.42  x  I0”6 

2.59  x  IO”4 

3.42  x  IO”3 

0.02 

0.04 

0 

0 

8.95  x  10"  11 

8.83  x  10”6 

7.71  x  10”  4 

7.12  x  10”  3 

0.03 

0.06 

0 

0 

2.88  x  IO-9 

4.79  x  10”  5 

2.12  x  10“  3 

0.01 

0.04 

0.09 

0 

0 

6.95  x  10“  8 

221  x  10”4 

5.42  x  10”  3 

0.03 

0.07 

0.12 

0 

0 

1.26  x  10-6 

9.45  x  10“4 

0.01 

0.05 

0.1 

0.17 

0 

3.8  X  10” 

15 

1.73  x  IO-5 

3.44  x  10”  3 

0.03 

0.08 

0.16 

0.24 

0 

6.66  x  10 

12 

1.78  x  lo"4 

0.01 

0.06 

0.14 

0.23 

0.32 

v° 

4.34  x  10' 

-9 

1.39  x  10”  3 

0.03 

0.11 

0.21 

0.32 

0.42  , 

Apply  binary  filter 


'o 

0 

0 

0 

0 

0 

0 

0^ 

gans.  .  :=  G.  .  Final. 
6  J.i  J,i  J 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-  14 

_  .  7 

0 

0 

3.64  x  10 

2x  10 

0 

0 

0 

0 

-  12 

-6 

-  4 

-3 

0 

0 

2.09  x  10 

1.42  x  10 

2.59  x  10 

3.42  x  10 

0 

0 

0 

0 

8.95  x  10”  11 

8.83  x  10”6 

7.71  x  10“  4 

7.12  x  IO-3 

0.03 

0 

_  9 

-5 

-3 

0 

0 

2.88  x  10 

4.79  x  10 

2.12  x  10 

0.01 

0 

0 

-  8 

-4 

Maximum  integral 

0 

0 

6.95  x  10 

2.27  x  10 

5.42  x  10 

0 

0 

0 

-  6 

-4 

max  (gans)  =0.03 

0 

0 

1.26  x  10 

9.45  x  10 

0 

0 

0 

0 

0 

3.8  x  10”  15 

1.73  x  10”  5 

0 

0 

0 

0 

0 

0 

6.66  x  10”  12 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0. 

C-10 


10  sieve 


H  = 


hans  = 


/ 

-  11 

-  6 

-4 

-  \ 

> 

0 

0 

1.28  x  10 

3.33  x  10 

4.07  x  10 

4.33  x  10 

0.02 

0.04 

0 

0 

4.81  x  10“  10 

1.92  x  10"  5 

1.16  x  10"3 

8.72  x 

10"3 

0.03 

0.05 

-  8 

-  5 

,  -  3 

0 

0 

1.36  x  10 

9.73  x  10 

3.06  x  10 

0.02 

0.04 

0.08 

-7 

-  4 

-3 

0 

0 

2.86  x  10 

4.3  x  10 

7.47  x  10 

0.03 

0.07 

0.1 

-6 

,  3 

0 

0 

4.51  x  10 

1.67  x  10 

0.02 

0.05 

0.1 

0.14 

-  13 

-5 

-  3 

0 

1.18  X  10 

5.36  x  10 

5.66  x  10 

0.04 

0.09 

0.14 

0.18 

-  10 

-4 

0 

1.38  x  10 

4.79  x  10 

0.02 

0.07 

0.14 

0.2 

0.24 

-  8 

-  3 

0 

6.05  x  10 

3.24  x  10 

0.04 

0.13 

0.2 

0.26 

0.29 

-  j 

0 

1.01  x  10 

0.02 

0.1 

0.21 

0.29 

0.34 

0.36 

-  4 

0 

6.42  x  10 

0.07 

0.21 

0.33 

0.39 

0.42 

0.42 

0 

0.02 

0.2 

0.38 

0.48 

0.51 

0.5 

0.48 

,0 

0.16 

0.46 

0.61 

0.65 

0.63 

0.58 

0.54  , 

Apply  binary  filter 

' 0 

0 

0 

0 

0  0 

0 

0^ 

nans. 

:■*  H  Final 

0 

0 

0 

0 

0  0 

0 

0 

i 

1  j.'  j.' 

-  8 

0 

0 

1.36  x  10 

0 

0  0 

0 

0 

-7 

-  4 

0 

0 

2.86  x  10 

4.3  x  10 

0  0 

0 

0 

-  6 

_  3 

0 

0 

4.51  x  10 

1.67  x  10 

0.02  0.05 

0 

0 

-  13 

-  5 

-3 

0 

1.18  x  10 

5.36  x  10 

5.66  x  10 

0.04  0.09  0.14 

0 

-  10 

-  4 

Maximum  integral 

0 

1.38  x  10 

4.79  x  10 

0.02 

0.07  0.14 

0 

0 

max(hans)  =  0.14 

„  _  -  8 

-  3 

0 

6.05  x  10 

3.24  x  10 

0.04 

0.13  0 

0 

0 

-  5 

0 

1.01  x  10 

0.02 

0.1 

0  0 

0 

0 

-4 

0 

6.42  x  10 

0.07 

0 

0  0 

0 

0 

0 

0.02 

0 

0 

0  0 

0 

0 

v0 

0 

0 

0 

0  0 

0 

0 , 

C-11 


14  sieve 


'  1 .41 

0.97 

0.88 

0.85 

1.46 

1.36 

1.17 

1.08 

1.51 

1.5 

1.39 

1.27 

1.55 

1.55 

1.51 

1.43 

1.6 

1.6 

1.59 

1.54 

1.65 

1.65 

1.65 

1.62 

1.7 

1.7 

1.7 

1.68 

1.75 

1.75 

1.75 

1.71 

1.8 

1.8 

1.79 

1.71 

1.86 

1.86 

1.8 

1.65 

1.92 

1.9 

1.72 

1.53 

^  1 .97 

1.81 

1.51 

1.34 

0.84 

0.84 

0.84 

0.82' 

1.03 

0.99 

0.96 

0.92 

1.19 

1.13 

1.07 

1.01 

1.34 

1.26 

1.17 

1.09 

1.46 

1.36 

1.26 

1.16 

1.54 

1.43 

1.31 

\2 

1.59 

1.47 

1.34 

1.22 

1.61 

1.48 

05 

1.22 

1.58 

1.45 

1.32 

\2 

1.51 

1.38 

1.26 

1.16 

1.39 

1.28 

1.18 

1.09 

1.23 

1.15 

1.08 

1.01, 

'1.41 

0 

0 

0 

0 

0 

0 

0^ 

1.46 

1.36 

0 

0 

0 

0 

0 

0 

1.51 

1.5 

1.39 

0 

0 

0 

0 

0 

1.55 

1.55 

1.51 

1.43 

0 

0 

0 

0 

1.6 

1.6 

1.59 

1.54 

1.46 

1.36 

0 

0 

1.65 

1.65 

1.65 

1.62 

1.54 

1.43 

1.31 

0 

1.7 

1.7 

1.7 

1.68 

1.59 

1.47 

0 

0 

1.75 

1.75 

1.75 

1.71 

1.61 

0 

0 

0 

1.8 

1.8 

1.79 

1.71 

0 

0 

0 

0 

1.86 

1.86 

1.8 

0 

0 

0 

0 

0 

1.92 

1.9 

0 

0 

0 

0 

0 

0 

J.97 

0 

0 

0 

0 

0 

0 

0, 

Apply  binary  filter 

kans.  .  :=  K.  .  Final. 
J.i  J.i  J> 


Maximum  integral 
max(kans)  =  1.97 


C-12 


30  sieve 


L 


Ians  = 


4.77  x  10”3 

0.45 

0.54 

0.57 

0.58 

0.58 

0.58 

0.58 

0 

0.1 

0.29 

0.39 

0.44 

0.47 

0.49 

0.5 

0 

9.88  x  10“  3 

0.12 

0.24 

0.32 

0.37 

0.4 

0.43 

0 

3.98  x  10"4 

0.04 

0.13 

0.21 

0.28 

0.32 

0.36 

0 

6.46  x  10"6 

0.01 

0.06 

0.14 

0.2 

0.25 

0.29 

0 

4.12  x  10"  8 

2.03  x  10" 

3 

0.03 

0.08 

0.14 

0.19 

0.23 

0 

1.01  x  10"  10 

-  4 

3  x  10 

0.01 

0.04 

0.09 

0.14 

0.18 

0 

-  14 

9.32  x  10 

3.34  x  10" 

5 

3.56  x  10"  3 

0.02 

0.06 

0.1 

0.14 

0 

0 

2.81  x  10" 

6 

1.05  x  10"  3 

0.01 

0.03 

0.07 

0.1 

0 

0 

1.77  x  10" 

7 

2.7  x  10"4 

4.7  x  10"  3 

0.02 

0.04 

0.07 

0 

0 

8.43  x  10" 

9 

6.11  x  10”  5 

1.92  x  10"  3 

0.01 

0.03 

0.05 

0 

0 

,  ,„-io 

3x10 

1.21  x  10"  5 

7.28  x  10”  4 

5.54  x  10"  3 

0.02 

0.04 

x  10  3 

0 

0 

0 

0 

0 

0 

0 

0 

0.1 

0 

0 

0 

0 

0 

0 

0 

9.88  x  10“  3 

0.12 

0 

0 

0 

0 

0 

0 

3.98  x  10"  4 

0.04 

0.13 

0 

0 

0 

0 

0 

6.46  x  10”  6 

0.01 

0.06 

0.14 

0.2 

0 

0 

0 

4.12  x  10"  8 

2.03  x  10"  3 

0.03 

0.08 

0.14 

0.19 

0 

0 

1.01  x  10"  10 

-  4 

3  x  10 

0.01 

0.04 

0.09 

0 

0 

0 

-  14 

9.32  x  10 

3.34  x  10"  5 

3.56  x  10"  3 

0.02 

0 

0 

0 

0 

0 

2.81  x  10"6 

1.05  x  10"  3 

0 

0 

0 

0 

0 

0 

1.77  x  10"7 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0, 

Apply  binary  filter 
Ians.  .  :=  L.  .  Final. 

I.i  I'  J 


Maximum  integral 
max(lans)  =  0.2 


V 


C-13 


Examples  of  expected  retained  fractions  across  NEDU  mesh  sizes  for  distributions  meeting 
the  overall  specifications.  The  output  file  of  this  MathCad  program  contains  46  such  distributions.  The 
Visual  Basic  program  "MeshFit"  fits  NEDU's  data  from  the  sieve  analysis  to  these  distributions.  The 
model  distribution  representing  the  best  fit  to  the  test  sample  is  then  compared  to  the  sample 
distribution.  That  comparison  is  quantified  by  assessments  of  "Goodness  of  Fit"  using  either  a  simple 
summed  error,  or  the  Chi  Square  statistical  test  (see  MeshFit  documentation). 


Definition  of  one  of  the  46  model  distributions  -  down  selected  from  135  candidate 
distributions. 


j:=3 

i  :=  1 


mmMean  :=  1.3  +  (0.05-j) 
mmMean  =  1.45 
mmSD  :=  (.05  -i) 


MeshMean  :=  1.096  + 


mmMean 


mmSD  =  0.05 


McsbMean  =  12.75 


8  mesh 


Fractions  retained  in  the  NEDU  sieves 
10  mesh  14  mesh  30  mesh 


ganSj  .  =  0 


hans.  .  =  0 
J.' 


kans.  .  s  1.51 

j.i 


Ians.  .  =  0 
J.' 


Another  model  distribution 
j  :=  12  mmMean  :=  1.3  +  (0.05-j)  MeshMean  :=  1.096+  I  16.9 


1 


mmMean 


MeshMean  =  9.99 


mmMean  =  1.9 
i  :=  1  mmSD  :=  (.05 -i) 

8  mesh  10  mesh 


mmSD  =  0.05 


ganSj  i  =  ° 


hans.  .  =  0 
J.i 


14  mesh 

kans.  .  =  1.97 
J.i 


30  mesh 

Ians.  ,  =  0 
J.I 


Yet  another  model  distribution 
j  :=  9  mmMean  :=  1.3  +  (0.05 •  j)  MeshMean  :=  1.096  +  I  16.9' 


1 


mmMean 


mmMean  =  1.75 
i  :=  8  mmSD  :=  (.05-i) 


mmSD  =  0.4 


MeshMean  =  10.75 


ganSj  .  =  0 


hans.  .  =  0 
J.i 


kanSj  ,  =  0 


Ians.  .  =  0 
J.' 


j  :=  6  mmMean  :=  1.3  +  (0.05-j)  MeshMean  :=  1.096  +  I  16.9 


mmMean  =  1.6 
i  :=  4  mmSD  :=  (.05-i) 


^16.9 


1 


mmMean  J 


minSD  =  0.2 


ganSj  .  =  8.83  x  10  6  hanSj  j  =  5.66  x  10 


kans.  =  1 .62 

j.i 


MeshMean  =  11.66 


Ians,  ,  =  0.03 
J  »  ' 


C-14 


j:~  12  mmMcan  :=  1.3+(0.05j)  McshMean  :=  1.096+  16.9- 


mmMean  =  1.9 
i  :=  4  mmSD  :=  (.05 -i) 


1 


mmMean  J 


McshMean  =  9.99 


mmSD  =  0.2 


ganSj  ,  =  0 


hans.  .  =  0 
J.' 


kans.  .  =  0 

J,i 


Ians.  .  =  0 
J.' 


j:=9  mmMean  :=  1.3  +  (0.05 -j)  MeshMean  :=  1.096+  16.9 


mmMean  =  1.75 
i  :=  4  mmSD  :=  (,05-i) 


1 


mmMean  J 


MeshMean  =  10.75 


mmSD  =  0.2 


-  4  -  3 

gans.  .  =  9.45  x  10  hans.  .  =  0.1  kans.  .  =  1.71  Ians.  .  =  1.05  x  10 

j,i  J,  •  J.'  J.' 


C-15 


APPENDIX  D.  Details  listing  for  Figure  10. 


no . 

Error 

3 . 5  mesh 

5  mesh 

8  mesh 

10  mesh 

14  mesh 

30  mesh 

i 

.  995 

0 

0 

0 

0 

0 

0 

2 

.  995 

0 

0 

0 

0 

0 

0 

3 

.  995 

0 

0 

0 

0 

0 

0 

4 

.  995 

0 

0 

0 

0 

0 

0 

5 

.  995 

0 

0 

0 

0 

0 

0 

6 

.  995 

0 

0 

0 

0 

0 

0 

7 

.  995 

0 

0 

0 

0 

0 

0 

8 

.  995 

0 

0 

0 

0 

0 

0 

9 

.3332 

5 . 771E-05 

0.07842 

0.7875 

0.0757 

0.05342 

0.004806 

10 

.  995 

0 

0 

0 

0 

0 

0 

11 

.  995 

0 

0 

0 

0 

0 

0 

12 

.  995 

0 

0 

0 

0 

0 

0 

13 

.  995 

0 

0 

0 

0 

0 

0 

14 

.  995 

0 

0 

0 

0 

0 

0 

15 

.  995 

0 

0 

0 

0 

0 

0 

16 

.  995 

0 

0 

0 

0 

0 

0 

17 

.30305 

8 . 588E-05 

0.09411 

0.7922 

0.06605 

0.04391 

0.003586 

18 

.26701 

0.0002359 

0 .1107 

0.758 

0.07048 

0.05423 

0.006197 

19 

.23123 

0.0005418 

0 . 1266 

0.7252 

0.07363 

0.064 

0.009646 

20 

.  995 

0 

0 

0 

0 

0 

0 

21 

.  995 

0 

0 

0 

0 

0 

0 

22 

.  995 

0 

0 

0 

0 

0 

0 

23 

.  995 

0 

0 

0 

0 

0 

0 

24 

.  995 

0 

0 

0 

0 

0 

0 

25 

.26821 

0.0001266 

0 . 112 

0.7924 

0.05706 

0.03575 

0.00265 

26 

.23136 

0.0003313 

0 . 1292 

0.7581 

0.06211 

0.04535 

0.004763 

27 

.19563 

0.0007314 

0 . 1454 

0.7254 

0.06594 

0.05469 

0.007653 

28 

.16117 

0.001406 

0 .1603 

0 . 6943 

0.06872 

0.06345 

0.01128 

29 

.  995 

0 

0 

0 

0 

0 

0 

30 

.  995 

0 

0 

0 

0 

0 

0 

31 

.  995 

0 

0 

0 

0 

0 

0 

32 

.  995 

0 

0 

0 

0 

0 

0 

33 

.23082 

0.0001851 

0 . 1322 

0.788 

0.04882 

0.02884 

0.00194 

34 

.19148 

0.0004617 

0 . 1497 

0.7543 

0.05427 

0.03761 

0.003631 

35 

.  15603 

0.0009675 

0 .1659 

0.7219 

0.05862 

0.04641 

0.006029 

36 

.12242 

0.001829 

0 .1806 

0 . 6912 

0.06195 

0.05487 

0.009132 

37 

.13617 

0.003116 

0.1937 

0 . 6623 

0.06442 

0.06275 

0.01287 

38 

.  995 

0 

0 

0 

0 

0 

0 

D-l 


39  .995  000000 

40  .995  000000 


41 

no . 

.20119 

Error 

0.000268 

3 . 5  mesh 

0 . 1547 

5  mesh 

0.7792 

8  mesh 

0.04135 

10  mesh 

0.02304 

14  mesh 

0.001407 

30  mesh 

42 

.  15302 

0.0006385 

0 . 1723 

0.7464 

0.04702 

0.03094 

0.002745 

43 

.11262 

0.001292 

0 . 1881 

0.7149 

0.05172 

0.0391 

0.004716 

44 

.08034 

0.002364 

0.2023 

0.685 

0.05549 

0.04715 

0.007347 

45 

. 10663 

0.00392 

0.2148 

0.6568 

0.05842 

0.05481 

0.0106 

46 

.13164 

0.006017 

0.2257 

0.6302 

0.06062 

0.06192 

0.01439 

47 

. 15536 

0.008681 

0.2349 

0.6053 

0.0622 

0.06836 

0.01861 

48 

.  995 

0 

0 

0 

0 

0 

0 

49 

.16505 

0.0003845 

0 . 1796 

0.7661 

0.03469 

0.01824 

0.001011 

50 

.  11765 

0.0008571 

0 .1968 

0.7346 

0.0404 

0.02525 

0.002059 

51 

.0747 

0.001713 

0.2121 

0.7044 

0.0453 

0.0327 

0.003663 

52 

.05458 

0.003037 

0.2255 

0.6758 

0.04938 

0.04026 

0.005874 

53 

.07959 

0.004906 

0.2371 

0 . 6486 

0.05267 

0.04762 

0.008683 

54 

. 10885 

0.007366 

0.2469 

0 . 6229 

0.05527 

0.05458 

0.01203 

55 

.  995 

0 

0 

0 

0 

0 

0 

56 

.  995 

0 

0 

0 

0 

0 

0 

57 

.12368 

0.0005469 

0.2069 

0.7487 

0.02881 

0.01431 

0.0007197 

58 

.07683 

0.001172 

0.2233 

0.7191 

0.03441 

0.02044 

0.001532 

59 

.03588 

0.002255 

0.2376 

0.6907 

0.03938 

0.02716 

0.002825 

60 

.05191 

0.003878 

0.2499 

0.6635 

0.04365 

0.03416 

0.004667 

61 

.09162 

0.006107 

0.2604 

0.6377 

0.04722 

0.04114 

0.007073 

62 

.  995 

0 

0 

0 

0 

0 

0 

63 

.  995 

0 

0 

0 

0 

0 

0 

64 

.  995 

0 

0 

0 

0 

0 

0 

65 

.08128 

0.0007443 

0.2364 

0.7275 

0.0237 

0.01112 

0.0005074 

66 

.04414 

0.001589 

0.2517 

0.7001 

0.02907 

0.01641 

0.00113 

67 

.05431 

0.002949 

0.2647 

0 . 6738 

0.03399 

0.0224 

0.002163 

68 

.09256 

0.004922 

0.2756 

0 . 6485 

0.03834 

0.0288 

0.003685 

69 

.  995 

0 

0 

0 

0 

0 

0 

70 

.  995 

0 

0 

0 

0 

0 

0 

71 

.  995 

0 

0 

0 

0 

0 

0 

72 

.  995 

0 

0 

0 

0 

0 

0 

73 

.07784 

0.001046 

0.2681 

0.7026 

0.0193 

0.008556 

0.0003544 

74 

.07139 

0.002137 

0.2817 

0 . 6779 

0.02434 

0.01307 

0.0008274 

D-2 


75 

.  10333 

0.00383 

0.2931 

0.654 

0.02911 

0.01834 

0. 

.001645 

76 

.  995 

0 

0 

0 

0 

0 

0 

77 

.  995 

0 

0 

0 

0 

0 

0 

78 

.  995 

0 

0 

0 

0 

0 

0 

79 

.  995 

0 

0 

0 

0 

0 

0 

80 

.  995 

0 

0 

0 

0 

0 

0 

81 

0.0002452 

.  11011 

0.001456 

0.3017 

0.6745 

0.01557 

0.006524 

82 

.13507 

0.002852 

0.3133 

0.6527 

0.02021 

0.01032 

0.0006007 

0 

84 

83 

.  995  0 

0 

0 

0 

0 

.  995 

0 

0 

0 

0 

0 

0 

no . 

Error 

3 . 5  mesh 

5  mesh 

8  mesh 

10  mesh 

14  mesh 

30 

mesh 

85 

.  995 

0 

0 

0 

0 

0 

0 

86 

.  995 

0 

0 

0 

0 

0 

0 

87 

.  995 

0 

0 

0 

0 

0 

0 

88 

.  995 

0 

0 

0 

0 

0 

0 

89 

0.0001681 

.18104 

0.002007 

0.3369 

0 . 6435 

0.01243 

0.004927 

90 

.  995 

0 

0 

0 

0 

0 

0 

91 

.  995 

0 

0 

0 

0 

0 

0 

92 

.  995 

0 

0 

0 

0 

0 

0 

93 

.  995 

0 

0 

0 

0 

0 

0 

94 

.  995 

0 

0 

0 

0 

0 

0 

95 

.  995 

0 

0 

0 

0 

0 

0 

96 

.  995 

0 

0 

0 

0 

0 

0 

Best  model 

=  59 

3 . 5  mesh 

5  mesh 

8  mesh 

10  mesh 

14  mesh 

30  mesh 

obs : 

0 

24.5 

68.5 

7 

0 

0 

modi : 

0 

23.9855 

69.07 

6.9365 

0 

0 

Chi  Square 

=.0163 

degrees  of 

freedom  =  2 

D-3 


APPENDIX  E.  Details  listing  for  Figure  11. 


no . 

Error 

3 . 5  mesh 

5  mesh 

8  mesh 

10  mesh 

14  mesh 

30 

i 

.  9916 

0 

0 

0 

0 

0 

0 

2 

.  9916 

0 

0 

0 

0 

0 

0 

3 

.  9916 

0 

0 

0 

0 

0 

0 

4 

.  9916 

0 

0 

0 

0 

0 

0 

5 

.  9916 

0 

0 

0 

0 

0 

0 

6 

.  9916 

0 

0 

0 

0 

0 

0 

7 

.  9916 

0 

0 

0 

0 

0 

0 

8 

.  9916 

0 

0 

0 

0 

0 

0 

9 

.  9482 

5 . 771E-05 

0.07842 

0.7875 

0.0757 

0.05342 

0 

10 

.  9916 

0 

0 

0 

0 

0 

0 

11 

.  9916 

0 

0 

0 

0 

0 

0 

12 

.  9916 

0 

0 

0 

0 

0 

0 

13 

.  9916 

0 

0 

0 

0 

0 

0 

14 

.  9916 

0 

0 

0 

0 

0 

0 

15 

.  9916 

0 

0 

0 

0 

0 

0 

16 

.  9916 

0 

0 

0 

0 

0 

0 

17 

.  91805 

8 . 588E-05 

0.09411 

0.7922 

0.06605 

0.04391 

0 

18 

.88201 

0.0002359 

0 .1107 

0.758 

0.07048 

0.05423 

0 

19 

.  84623 

0.0005418 

0 . 1266 

0.7252 

0.07363 

0.064 

0 

20 

.  9916 

0 

0 

0 

0 

0 

0 

21 

.  9916 

0 

0 

0 

0 

0 

0 

22 

.  9916 

0 

0 

0 

0 

0 

0 

23 

.  9916 

0 

0 

0 

0 

0 

0 

24 

.  9916 

0 

0 

0 

0 

0 

0 

25 

.88321 

0.0001266 

0 . 112 

0.7924 

0.05706 

0.03575 

0 

26 

.  84636 

0.0003313 

0 . 1292 

0.7581 

0.06211 

0.04535 

0 

27 

.81063 

0.0007314 

0 . 1454 

0.7254 

0.06594 

0.05469 

0 

28 

.77617 

0.001406 

0 .1603 

0 . 6943 

0.06872 

0.06345 

0 

29 

.  9916 

0 

0 

0 

0 

0 

0 

30 

.  9916 

0 

0 

0 

0 

0 

0 

31 

.  9916 

0 

0 

0 

0 

0 

0 

32 

.  9916 

0 

0 

0 

0 

0 

0 

33 

.  84346 

0.0001851 

0 . 1322 

0.788 

0.04882 

0.02884 

0 

34 

.  80648 

0.0004617 

0 . 1497 

0.7543 

0.05427 

0.03761 

0 

35 

.77103 

0.0009675 

0 .1659 

0.7219 

0.05862 

0.04641 

0 

36 

.73742 

0.001829 

0 .1806 

0 . 6912 

0.06195 

0.05487 

0 

37 

.70577 

0.003116 

0 .1937 

0 . 6623 

0.06442 

0.06275 

0 

38 

.  9916 

0 

0 

0 

0 

0 

0 

39 

.  9916 

0 

0 

0 

0 

0 

0 

E-l 


mesh 


.004806 


.003586 
.006197 
.009646 


.00265 
.004763 
.007653 
.01128 


.00194 
.003631 
.006029 
.009132 
.01287 


40 

41 

42 

no . 

.  9916 
.79889 

.76206 

Error 

0 

0.000268 

0.0006385 

3 . 5  mesh 

0 

0 . 1547 

0 . 1723 

5  mesh 

0 

0.7792 

0.7464 

8  mesh 

0 

0.04135 

0.04702 

10  mesh 

0 

0.02304 

0.03094 

14  mesh 

0 

0 

0 

30 

43 

.72762 

0.001292 

0 .1881 

0.7149 

0.05172 

0.0391 

0 

44 

.  69534 

0.002364 

0.2023 

0.685 

0.05549 

0.04715 

0 

45 

.  66523 

0.00392 

0.2148 

0.6568 

0.05842 

0.05481 

0 

46 

.  63704 

0.006017 

0.2257 

0.6302 

0.06062 

0.06192 

0 

47 

.  61096 

0.008681 

0.2349 

0.6053 

0.0622 

0.06836 

0 

48 

.  9916 

0 

0 

0 

0 

0 

0 

49 

.74943 

0.0003845 

0 . 1796 

0.7661 

0.03469 

0.01824 

0 

50 

.71345 

0.0008571 

0 .1968 

0.7346 

0.0404 

0.02525 

0 

51 

.  6803 

0.001713 

0.2121 

0.7044 

0.0453 

0.0327 

0 

52 

.  64994 

0.003037 

0.2255 

0.6758 

0.04938 

0.04026 

0 

53 

.  62179 

0.004906 

0.2371 

0.6486 

0.05267 

0.04762 

0 

54 

.59585 

0.007366 

0.2469 

0.6229 

0.05527 

0.05458 

0 

55 

.  9916 

0 

0 

0 

0 

0 

0 

56 

.  9916 

0 

0 

0 

0 

0 

0 

57 

.  69492 

0.0005469 

0.2069 

0.7487 

0.02881 

0.01431 

0.0007197 

58 

.  66065 

0.001172 

0.2233 

0.7191 

0.03441 

0.02044 

0 

59 

.  62964 

0.002255 

0.2376 

0.6907 

0.03938 

0.02716 

0 

60 

.  60141 

0.003878 

0.2499 

0.6635 

0.04365 

0.03416 

0 

61 

.57566 

0.006107 

0.2604 

0.6377 

0.04722 

0.04114 

0 

62 

.  9916 

0 

0 

0 

0 

0 

0 

63 

.  9916 

0 

0 

0 

0 

0 

0 

64 

.  9916 

0 

0 

0 

0 

0 

0 

65 

.  63592 

0.0007443 

0.2364 

0.7275 

0.0237 

0.01112 

0.0005074 

66 

.  60388 

0.001589 

0.2517 

0.7001 

0.02907 

0.01641 

0 

67 

.57549 

0.002949 

0.2647 

0.6738 

0.03399 

0.0224 

0 

68 

.55004 

0.004922 

0.2756 

0.6485 

0.03834 

0.0288 

0 

69 

.  9916 

0 

0 

0 

0 

0 

0 

70 

.  9916 

0 

0 

0 

0 

0 

0 

71 

.  9916 

0 

0 

0 

0 

0 

0 

72 

.  9916 

0 

0 

0 

0 

0 

0 

73 

.57236 

0.001046 

0.2681 

0.7026 

0.0193 

0.008556 

0.0003544 

74 

.54361 

0.002137 

0.2817 

0.6779 

0.02434 

0.01307 

0.0008274 

75 

.51835 

0.00383 

0.2931 

0.654 

0.02911 

0.01834 

0 

E-2 


.001407 

.002745 

mesh 

.004716 

.007347 

.0106 

.01439 

.01861 

.001011 
.002059 
.003663 
.005874 
.008683 
.01203 


.001532 
.002825 
.004667 
.007073 


.00113 
.002163 
.003685 


.001645 


76 

.  9916 

0 

0 

0 

0 

0 

0 

77 

.  9916 

0 

0 

0 

0 

0 

0 

78 

.  9916 

0 

0 

0 

0 

0 

0 

79 

.  9916 

0 

0 

0 

0 

0 

0 

80 

.  9916 

0 

0 

0 

0 

0 

0 

81 

0.0002452 

.50489 

0.001456 

0.3017 

0.6745 

0.01557 

0.006524 

82 

0.0006007 

.47993 

0.002852 

0.3133 

0.6527 

0.02021 

0.01032 

83 

.  9916 

0 

0 

0 

0 

0 

0 

84 

.  9916 

0 

0 

0 

0 

0 

0 

85 

.  9916 

0 

0 

0 

0 

0 

0 

no . 

Error 

3 . 5  mesh 

5  mesh 

8  mesh 

10  mesh 

14  mesh 

30  mesh 

86 

.  9916 

0 

0 

0 

0 

0 

0 

87 

.  9916 

0 

0 

0 

0 

0 

0 

88 

.  9916 

0 

0 

0 

0 

0 

0 

89 

.4367 

0.002007 

0.3369 

0.6435 

0.01243 

0.004927 

0.0001681 

90 

.  9916 

0 

0 

0 

0 

0 

0 

91 

.  9916 

0 

0 

0 

0 

0 

0 

92 

.  9916 

0 

0 

0 

0 

0 

0 

93 

.  9916 

0 

0 

0 

0 

0 

0 

94 

.  9916 

0 

0 

0 

0 

0 

0 

95 

.  9916 

0 

0 

0 

0 

0 

0 

96 

.  9916 

0 

0 

0 

0 

0 

0 

Best  model 

3 . 5  mesh 

=  89 

5  mesh 

8  mesh 

10  mesh 

14  mesh 

30  mesh 

obs :  0 

55.1 

44.12 

0 

0 

0 

modi :  0 

33.8907 

66.10251 

0 

0 

0 

Chi  Square 

=  20.5834 

degrees  of 

freedom  =  1 

E-3 


APPENDIX  F.  MeshFit  Code 

VERSION  4.00 

Dim  f racRetain ( 8 ,  50,  50)  As  Single 

Dim  Frac3  As  Single 

Dim  Frac5  As  Single 

Dim  Frac8  As  Single 

Dim  FraclO  As  Single 

Dim  Fracl4  As  Single 

Dim  Frac30  As  Single 

Dim  mesh35fit  As  Single 

Dim  mesh5fit  As  Single 

Dim  mesh8fit  As  Single 

Dim  meshlOfit  As  Single 

Dim  meshl4fit  As  Single 

Dim  mesh30fit  As  Single 

Dim  ErrorFrac  As  Single 
Dim  Fractotal  As  Single 

Dim  i  As  Integer 

Dim  imax  As  Integer 

Dim  j  As  Integer 

Dim  jmax  As  Integer 

Dim  s  As  Integer 

Dim  sieveNumb  As  Integer 

Dim  bestn  As  Integer 

Private  Sub  Commandl_Click ( ) 

If  Fractotal  >  1#  Then 
Beep 
Beep 
Beep 

Label ( 15 ). Visible  =  True 
Exit  Sub 
End  If 

Call  TotalFraction 

Description$  =  Text7.Text 

If  Choice . Optionl . Value  =  True  Then 

Open  "e:\meshfit\mesh408.prn"  For  Input  As  #1 
jmax  =  12 
imax  =  8 
End  If 

If  Choice . Option2 . Value  =  True  Then 

Open  "e:\meshfit\meshl014.prn"  For  Input  As  #1 
jmax  =  11 
imax  =  9 
End  If 

If  Choice . Option3 . Value  =  True  Then 

Open  "e:\meshfit\mesLn408.prn"  For  Input  As  #1 
jmax  =  12 
imax  =  8 
End  If 

j  =  1 

i  =  1 
s  =  1 

sieveNumb  =  6 
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Do  Until  s  >  sieveNumb 

j  =  1 

Do  Until  j  >  jmax 
i  =  1 

Do  Until  i  >  imax 

Input  #1,  fracRetain (s,  j,  i)  '  Get  data 
i  =  i  +  1 
Loop 

j  =  j  +  1 

Loop 

s  =  s  t  1 
Loop 

Close  #1  '  Close  file, 

minerror  =  50 

n  =  1 

If  Checkl . Value  =  1  Then 

Printer . Print  'create  blank  line 

Printer . Font . Name  =  "Arial" 

Printer . Font . Size  =  10 

Printer . Print  Description$  'print  test  description 
Printer . Print 

Printer . Font . Underline  =  True 

Printer . Print  "no.",  "Error",  "3.5  mesh",  "5  mesh",  "8  mesh",  "10  mesh",  "14 
mesh",  "30  mesh" 

Printer . Font . Under line  =  False 
End  If 


For 


j  = 

For 


1  To  jmax 
i  =  1  To  imax 

If  Choice . Optionl . Value  =  True 
distMean  =  3.1  t  (0.05  *  j) 
distSD  =  0.6  +  (0.05  *  i) 

Else 

distMean  =  1.3  t  (0.05  *  j) 
distSD  =  0.05  *  i 
End  If 


Then  ' for  408 

'mean  in  mm 
'SD  in  mm 
' for  812  mesh 
'mean  in  mm 
'SD  in  mm 


mesh 


ErrorFrac35  =  Abs(Frac35 
ErrorFrac5  =  Abs (Frac5  - 
ErrorFrac8  =  Abs (Frac8  - 
ErrorFraclO  =  Abs(Fracl0 
ErrorFracl4  =  Abs(Fracl4 
ErrorFrac30  =  Abs(Frac30 


-  fracRetain (1 ,  j,  i)) 
fracRetain (2,  j,  i)) 
fracRetain (3,  j,  i)) 

-  fracRetain (4 ,  j,  i)) 

-  fracRetain  (5,  j,  i)) 

-  fracRetain  ( 6 ,  j,  i)) 


sumerror  =  ErrorFrac5  t  ErrorFrac8  t  ErrorFraclO  t  ErrorFracl4 


sumerrorp  =  Format ( sumerror ,  "#.#####") 

If  Checkl. Value  =  1  Then  Printer . Print  n,  sumerrorp,  fracRetain ( 1 , 

j,  i),  fracRetain (2,  j,  i) ,  fracRetain (3,  j,  i) ,  fracRetain ( 4 ,  j,  i),  fracRetain ( 5,  j, 
i) ,  fracRetain  (6,  j,  i) 

End  If 


If  sumerror  <  minerror  Then 
bestn  =  n 
bestj  =  j 
besti  =  i 

minerror  =  sumerror 
bestdistMean  =  distMean 
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mesh  size 
mesh  size 


bestdistSD  =  distSD 

meshMean  =  1.096  +  (16.91  *  1  /  bestdistMean) 

meshSDup  =  1.096  +  (16.91  *  1  /  (bestdistMean 

meshSDdn  =  1.096  +  (16.91  *  1  /  (bestdistMean 


meshMean  =  Format (meshMean,  "##.#") 
meshSDup  =  Format (meshSDup,  "##.#") 
meshSDdn  =  Format (meshSDdn,  "##.#") 
End  If 

n  =  n  +  1  'increment  counter 

Next  i 
Next  j 

Label2 . Caption  =  Format$ (bestdistSD,  "#.#") 


'mean  in  mesh  size 
-  bestdistSD) )  'mean 

+  bestdistSD) )  'mean 


+ 


SD  in 
SD  in 


mesh35fit  =  f racRetain ( 1 ,  bestj 
mesh5fit  =  f racRetain ( 2 ,  bestj, 
mesh8fit  =  fracRetain  (3,  bestj, 
meshlOfit  =  fracRetain ( 4 ,  bestj 
meshl4fit  =  fracRetain ( 5 ,  bestj 
mesh30fit  =  fracRetain ( 6 ,  bestj 


,  besti) 
besti) 
besti) 

,  besti) 
,  besti) 
,  besti) 


Label  (20)  . Caption  =  Format$ (fracRetain  (1, 
Label (21 ). Caption  =  Format$ ( fracRetain (2 , 
Label ( 22 ). Caption  =  Format$ (fracRetain (3, 
Label ( 23 ). Caption  =  Format$ (fracRetain (4, 
Label (24 ). Caption  =  Format$ (fracRetain (5, 
Label (25) . Caption  =  Format$ (fracRetain (6, 
For  idx  =  20  To  25 


bestj , 
bestj , 
bestj , 
bestj , 
bestj , 
bestj , 


besti) 

100, 

"#.###" 

besti) 

* 

100, 

"#.###" 

besti) 

•k 

100, 

"#.###" 

besti) 

■k 

100, 

"#.###" 

besti) 

■k 

100, 

"#.###" 

besti) 

-k 

100, 

"#.###" 

If  Label ( idx) . Caption  =  "."  Then  Label ( idx) . Caption  =  "0" 
Next  idx 


Call  ChiSquare 
Call  BarGraph 


Label (37) . Caption 
Label (39) . Caption 
Label (11)  . Caption 
Label ( 18 ) . Caption 
Label (40) .Caption 
End  Sub 


Format$ (minerror,  "#.#####") 
Str$ (bestn) 

Str$ (meshMean) 

Str$ (meshSDup) 

Str$ (meshSDdn) 


Private  Sub  Command2_Click ( ) 
End 

End  Sub 


Private  Sub  Command3_Click ( ) 

Mesh . BackColor  =  &HFFFFFF 
Checkl . BackColor  =  &HFFFFFF 
For  idx  =  0  To  39 

Label ( idx) . BackColor  =  &HFFFFFF  'control  array 
Next  idx 

Printer . DrawMode 
Printer . ScaleLeft 
Printer . ScaleTop 
PrintForm 
Printer . EndDoc 

Mesh . BackColor  =  &H8000000F 
Checkl . BackColor  =  &H8000000F 


=  13 

=  -(( Printer . Width  -  Mesh. Width)  /  2) 

=  -( (Printer . Height  -  Mesh. Height)  /  2) 
'print  graphic  screen 
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'control  array 


For  idx  =  0  To  39 

Label ( idx) . BackColor  =  &H8000000F 
Next  idx 
End  Sub 

Private  Sub  Command4_Click ( ) 

Commandl . SetFocus 
Choice. Show  0 
End  Sub 

Private  Sub  Command5_Click ( ) 

Textl.Text  =  "" 

Text2 . Text  =  "" 

Text3 . Text  =  "" 

Text4 . Text  =  "" 

Text5.Text  =  "" 

Text 6. Text  =  "" 

End  Sub 

Private  Sub  Form_Load() 

Label ( 1 6 ). Visible  =  False 
Graphl . RandomData  =  0 
Graphl . NumPoints  =  6 

Graphl . LeftTitle  =  "Fraction  Retained" 
Graphl . BottomTitle  =  "Mesh  Size" 

Graphl . NumSets  =  2 
End  Sub 


Private  Sub  VScrolll_Change ( ) 

Text7.Text  =  Format$ (VScrolll .Value, 
End  Sub 


Private  Sub  Optionl_Click ( ) 
End  Sub 


Private  Sub  0ption2_Click ( ) 
End  Sub 


Private  Sub  Textl_Change ( ) 

Label ( 1 6 ). Visible  =  False 
Call  TotalFraction 
End  Sub 


Private  Sub  Text2_Change ( ) 

Label ( 1 6 ). Visible  =  False 
Call  TotalFraction 
End  Sub 


Private  Sub  Text3_Change ( ) 

Label ( 1 6 ). Visible  =  False 
Call  TotalFraction 
End  Sub 


"##") 
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Private  Sub  Text4_Change ( ) 

Label ( 1 6 ). Visible  =  False 
Call  TotalFraction 
End  Sub 


Private  Sub  Text5_Change ( ) 

Label ( 1 6 ). Visible  =  False 
Call  TotalFraction 
End  Sub 


Private  Sub  Text6_Change ( ) 

Label ( 1 6 ). Visible  =  False 
Call  TotalFraction 
End  Sub 


Public  Sub  TotalFraction  ( ) 

Frac3  =  Val (Textl . Text)  /  100 
Frac5  =  Val (Text2. Text)  /  100 
Frac8  =  Val (Text3 . Text)  /  100 
FraclO  =  Val  (Text4 . Text)  /  100 
Fracl4  =  Val (Text5 . Text)  /  100 
Frac30  =  Val (Text6 . Text)  /  100 

Fractotal  =  Frac3  +  Frac5  +  Frac8  +  FraclO  +  Fracl4 
If  Fractotal  <=  1  Then 

Gaugel. Value  =  100  *  Fractotal 
Else:  Gaugel. Value  =  0 
End  If 

Label  ( 15)  . Caption  -  Format$(100  *  Fractotal,  "#.##") 
End  Sub 


Public  Sub  BarGraph ( ) 
Graphl . DrawMode  = 
Graphl . PatternData 
Graphl . YAxisMax  = 
Graphl .Autolnc  =  1 
Graphl . LabelText  = 
Graphl . LabelText  = 
Graphl . LabelText  = 
Graphl . LabelText  = 
Graphl . LabelText  = 
Graphl . LabelText  = 
Graphl . GraphData  = 
Graphl . GraphData  = 
Graphl . GraphData  = 
Graphl . GraphData  = 
Graphl . GraphData  = 
Graphl . GraphData  = 
Graphl . GraphData  = 
Graphl . GraphData  = 
Graphl . GraphData  = 
Graphl . GraphData  = 
Graphl . GraphData  = 
Graphl . GraphData  = 
Graphl . LegendText 
Graphl . LegendText 
Graphl . DrawMode  = 
End  Sub 


0 

=  5 

1 

"  3.5" 

"  5" 

"  8" 

"  10" 

"  14" 

"  30" 

Frac3 
Frac5 
Frac8 
FraclO 
Fracl 4 
Frac30 
mesh35fit 
mesh5fit 
mesh8fit 
meshlOfit 
meshl4fit 
mesh30fit 
=  "measured" 
=  "best  fit" 
2 


Frac30 
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Public  Sub  ChiSquareO 
combined  =  0 
mult  =  100 

Frac3mult  =  mult  *  Frac3 
Frac5mult  =  mult  *  Frac5 
Frac8mult  =  mult  *  Frac8 
FraclOmult  =  mult  *  FraclO 
Fracl4mult  =  mult  *  Fracl4 
Frac30mult  =  mult  *  Frac30 

If  Frac30mult  <  5  Then 
combine30  =  1 

Fracl4mult  =  Fracl4mult  +  Frac30mult 
meshl4fit  =  meshl4fit  +  mesh30fit 
Frac30mult  =  0 
mesh30fit  =  0 
End  If 

If  Fracl4mult  <  5  Then 
combinel4  =  1 

FraclOmult  =  FraclOmult  +  Fracl4mult 
meshlOfit  =  meshlOfit  +  meshl4fit 
Fracl4mult  =  0 
meshl4fit  =  0 
End  If 

If  FraclOmult  <  5  Then 
combinelO  =  1 

Frac8mult  =  Frac8mult  +  FraclOmult 
mesh8fit  =  mesh8fit  +  meshlOfit 
FraclOmult  =  0 
meshlOfit  =  0 
End  If 

If  Frac8mult  <  5  Then 
combine8  =  1 

FraclOmult  =  FraclOmult  +  Frac8mult 
meshlOfit  =  meshlOfit  +  mesh8fit 
Frac8mult  =  0 
mesh8fit  =  0 
End  If 

If  Frac3mult  <  5  Then 
combine35  =  1 

Frac5mult  =  Frac5mult  +  Frac3mult 
mesh5fit  =  mesh5fit  +  mesh35fit 
Frac3mult  =  0 
mesh35fit  =  0 
End  If 

If  Frac5mult  <  5  Then 
combine5  =  1 

Frac8mult  =  Frac8mult  +  Frac5mult 
mesh8fit  =  mesh8fit  +  mesh5fit 
Frac5mult  =  0 
mesh5fit  =  0 
End  If 
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mesh35mult 
mesh5mult  = 
mesh8mult  = 
meshlOmult 
meshl4mult 
mesh30mult 


=  mult  *  mesh35fit 
mult  *  mesh5fit 
mult  *  mesh8fit 
=  mult  *  meshlOfit 
=  mult  *  meshl4fit 
=  mult  *  mesh30fit 


If  Checkl . Value  =  1  Then 
Printer . Print 

Printer . Print  "Best  model  =";  bestn 
Printer . Print 

Printer . Font . Under line  =  True 

Printer . Print  "  3.5  mesh",  "5  mesh",  "8  mesh",  "10  mesh",  "14  mesh",  "30  mesh" 

Printer . Font . Under line  =  False 

Printer . Print  "obs :  Frac3mult,  Frac5mult,  Frac8mult,  FraclOmult,  Fracl4mult, 

Frac30mult 

Printer . Print  "modi:  mesh35mult,  mesh5mult,  mesh8mult,  meshlOmult,  meshl4mult, 

mesh30mult 
End  If 


combined  =  combine35  +  combine5  +  combine8  +  combinelO  +  combinel4  +  combine30 
degfreedom  =  6  -  combined  -  1  '6  sieves  -  sieves  combined  -  1 

If  mesh35fit  >  0.00001  Then 

chi35  =  (Frac3mult  -  mesh35mult)  A  2  /  mesh35mult 
Else:  chi35  =  0 
End  If 

If  mesh5fit  >  0.00001  Then 

chi5  =  (Frac5mult  -  mesh5mult)  A  2  /  mesh5mult 
Else:  chi5  —  0 
End  If 

If  mesh8fit  >  0.00001  Then 

chi8  =  (Frac8mult  -  mesh8mult)  A  2  /  mesh8mult 
Else:  chi8  =  0 
End  If 

If  meshlOfit  >  0.00001  Then 

chilO  =  (FraclOmult  -  meshlOmult)  A  2  /  meshlOmult 
Else:  chilO  =  0 
End  If 

If  meshl4fit  >  0.00001  Then 

chil4  =  (Fracl4mult  -  meshl4mult)  A  2  /  meshl4mult 
Else:  chil4  =  0 
End  If 

If  mesh30fit  >  0.00001  Then 

chi30  =  (Frac30mult  -  mesh30mult)  A  2  /  mesh30mult 
Else:  chi30  =  0 
End  If 

ChiStatistic  =  chi35  +  chi5  +  chi8  +  chilO  +  chil4  +  chi30 
Label  (26)  . Caption  =  Format$ (ChiStatistic,  "###.####") 

Label ( 2 9 ). Caption  =  degfreedom 

If  Checkl. Value  =  1  Then 
Printer . Print 

Printer . Print  "Chi  Square  =" ;  Format$ (ChiStatistic,  "###.####") 

Printer . Print  "degrees  of  freedom  =";  degfreedom 
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End  If 


significance  =  0 

If  degfreedom  =  1  And  ChiStatistic  >  3.84146  Then  significance  =  1 

If  degfreedom  =  2  And  ChiStatistic  >  5.99147  Then  significance  =  1 

If  degfreedom  =  3  And  ChiStatistic  >  7.81473  Then  significance  =  1 

If  degfreedom  =  4  And  ChiStatistic  >  9.48773  Then  significance  =  1 

If  degfreedom  =  5  And  ChiStatistic  >  11.0705  Then  significance  =  1 

If  significance  =  1  Then 

Label (30) . ForeColor  =  &H0& 

Label ( 30 ). Caption  =  "Significantly  different  from  best  model  !" 
End  If 

If  significance  =  0  Then 

Label  (30)  . ForeColor  =  &H0 

Label (30) . Caption  =  "Same  as  model" 

End  If 

End  Sub 
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Appendix  G  MeshFit  SUMMARY 


When  NEDU  evaluates  sodalime  quality,  we  determine,  among  other  things,  whether  or 
not  the  absorbent  meets  specifications  for  granule  size  distributions.  However,  because 
of  disparities  between  mesh  sizes  called  out  in  specifications  and  sieve  sizes  used  in 
testing,  a  simple  comparison  is  not  always  possible.  MeshFit  software  makes  those 
comparisons  simple  for  the  operator,  and  gives  a  clear  answer  as  to  whether  a  sodalime 
sample  meets  specifications. 

The  software  relies  on  the  mathematical  translation  of  sodalime  specifications  into 
specific  targets  for  sieve-based  sodalime  testing.  First,  mathematical  distributions  of 
granule  sizes  are  determined  that  meet  the  sodalime  specification.  Those  distributions, 
called  model  distributions,  are  then  subjected  to  a  mathematical  sieving  operation.  After 
the  test  bench  operator  enters  the  results  of  the  actual  physical  sieving  operation, 
MeshFit  compares  those  results  to  the  “sieved”  model  distributions.  If  the  physical 
sample  matches  one  of  the  model  distributions,  then  the  sample  meets  the  granule  size 
specification. 

It  is  possible  that  a  test  sample  will  match  none  of  the  model  distributions  particularly 
well.  Fortunately,  the  degree  of  match  or  mismatch  can  be  measured  and  used  as  an 
evaluation  factor  for  the  test  sample.  MeshFit  finds  the  one  model  distribution  that  best 
agrees  with  the  test  sample  data.  A  statistical  test  then  determines  whether  or  not  the 
best  fit  is  adequate;  that  is,  whether  the  actual  sample  and  the  model  are  in  close 
enough  agreement.  The  goodness  of  fit  between  the  test  data  and  the  best  model 
distribution  is  found  through  the  classical  chi  square  statistical  test  for  goodness  of  fit. 

When  using  MeshFit,  the  user  selects  whether  the  test  sample  is  from  a  4-8  mesh 
sample  or  10-14  U.S.  mesh  (8-12  British  mesh)  sample.  The  user  then  enters  the  mean 
retained  percentages  for  all  sieve  screens  obtained  from  the  sieving  of  five  or  more 
sodalime  samples.  The  software  next  identifies  the  best  fit  to  a  model  distribution,  and 
plots  in  bar  graph  form  the  mesh  distribution  of  both  the  sample  and  the  best  model. 
Finally,  the  chi  square  goodness  of  fit  statistic  is  computed,  with  significance  being 
determined  from  the  chi  square  and  the  degrees  of  freedom  for  the  test  (Figures  G-1 
and  G-2). 

The  chi  square  test  enables  the  testing  of  observed  versus  expected  values  for  multiple 
data  groupings.  In  this  case,  we  compare  the  observed  fraction  retained  by  a  particular 
sieve  screen  against  the  expected  fraction  for  the  closest  normal  distribution  meeting 
the  specifications. 
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Figure  G-1.  MeshFit  printout  for  a  sample  of  812  Sofnolime  that  meets  the  British 
specification  for  the  granule  size  distribution. 
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Figure  G-2.  MeshFit  printout  for  a  sample  of  812  Sofnolime  that  does  not  meet  the 
British  specification  for  granule  size  distribution. 
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